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ABSTEACT 

The  external  nucleate  pool-boiling  heat-transfer  coeffi- 
cient of  a  horizontal  smooth  copper  tube  in  R-114-oil 
mixtures  (0  to  10  percent  oil)  was  measured  for  heat  fluxes 
from  1  to  100  kW/m2  at  two  different  saturation  temperatures 
(-2-2  °C  and  6.7  °C) .  A  copper-nickel  tube  coated  with  the 
Onion  Carbide  "High  Flux"  coating  was  similarly  tested.  The 
High  Flux  coating  was  found  to  improve  the  heat-transfer 
coefficient  by  at  least  a  factor  of  7  in  oil- free  R-114.  Oil 
resulted  in  about  a  20  percent  reduction  of  the  heat- 
transfer  coefficient  of  the  High  Flux  surface  at  heat  fluxes 
less  than  30  ktf/m2  and  up  to  an  80  percent  reduction  at  heat 
fluxes  above  30  ktf/m2  with  greater  than  6  percent  oil. 
Under  all  conditions,  the  High  Flux  coated  tube  outperformed 
the    smooth    copper   tube. 


TABLE  OP  CONTENTS 


I.  INTRODUCTION 10 

A.       BACKGP.ODND 10 

3.       THESIS    OBJECTIVES       12 

II.  REVIEW    OF    REFRIGERANT-OIL    MIXTURE    BEHAVIOR    ....  13 

A.  NUCLEATE    ECILING    OF    REFRIGERANT-OIL 

MIXTURES    FROM    SMOOTH    TUBES    13 

1.  Physical    Properties      ....13 

2.  Saturation   Temperature    20 

3.  Tube    Diameter 21 

4.  Surface   Condition      23 

5.  Hydrostatic  Effect 23 

B.  NUCLEATE    ECILING    OF    REFRIGERANT-OIL 

MIXTURES  FROM  ENHANCED  SURFACES   23 

1.  Fins  and  Surfaces  with  Reentrant 

Cavities 24 

2.  Porous  Coatings 27 

III.  DESCRIPTION    OF   EXPERIMENTAL    APPARATUS       30 

A.  OVERALL    APPARATUS       30 

B.  OIL    SAMPLING    APPARATUS 34 

C.  BOILING    TUBE    CONSTRUCTION 35 

1.  Short    Tubes 35 

2.  Long    Tube 38 

D.  DATA    ACQUISITION    AND    REDUCTION 40 

E.  COMPUTER-CCNTROLLED    VALVE       42 

IV.  EXPERIMENTAL    PROCEDURES       45 

A.  INSTALLATION    OF    TUBE    IN    APPARATUS 45 

B.  GENERAL    OPERATION       46 


C.   SAMPLING  OIL 50 

V.  SYSTEM  OPERATION  AND  PROBLEMS   54 

A.  CIRCUMFERENTIAL  WALI  TEMPERATURE  VARIATION  .  .  54 

B.  AXIAL  TEMPERATURE  VARIATION  OF  LONG  TUBES   .  .  57 

C.  COMPUTER-CONTROLLED  VALVE   6  1 

VI.  RESULTS  AND  DISCUSSION 64 

A.  BOILING  PERFORMANCE  OF  SMOOTH  TUBE 64 

B.  BOILING  PERFORMANCE  OF  HIGH  FLUX  SURFACE  ...  68 

C.  COMPARISON  OF  HIGH  FLUX  TO  SMOOTH  TUBE 

BOILING  PERFORMANCE   73 

D.  EFFECT  OF  SATURATION  TEMPERATURE  ON 

BOILING  PERFORMANCE   76 

E.  LOCAL  OIL  SAMPLING  EFFORTS  77 

VII.  CONCLUSIONS  AND  RECOMMENDATIONS   81 

A.  CONCLUSIONS 8  1 

B.  RECOMMENDATIONS 82 

APPENDIX  A:   THERMOCOUPLE  CALIBRATION  84 

APPENDIX  E:   DATA  REDUCTION  PROGRAM  88 

APPENDIX  C:   COMPUTER-CONTROLLED  VALVE  PROGRAM   ....   109 
APPENDIX  D:   EXAMPLES  OF  REPRESENTATIVE  DATA  RUNS  ...   112 

APPENDIX  E:   UNCERTAINTY  ANALYSIS  114 

LIST  OF  REFERENCES 119 

INITIAL  DISTRIBUTION  LIST  122 


LIST  OF  TABLES 

1.  Summary    of  Short  Tube   Assembly  Methods 38 

2.  HP-3497A    Channel   Allocations    42 

3.  Summary   of   Data  Runs 47 

4.  Configuration    cf    SYS-2A    I/O   Channels 109 

5.  Uncertainty   Analysis    Terms    Using    8 
Thermocouples      116 

6.  Uncertainty   Analysis    Terms    Using    Center   4 
Thermocouples      118 


LIST    OF    FIGURES 

1.  1         Schematic  and    Scanning    Electron   Micrograph 

(500x)    of   High    Flux  Surface 11 

2.1  Effect   of    Oil    ce   Boiling    Coefficient    of    R- 114 

(from   Ref.    6) 14 

2.2  Variation    of    Heat-Transfer   Coefficient  of 

R-114-0il    Mixtures    (from    Ref.     13)     15 

2.3  Surface-Tension  Variation   of   R-114-0il   Mixtures 

(from  Ref.    6) 17 

2.4  Oil   Concentration    Gradient    in    a   Bubble   in 
Refrigerant-Oil  Mixtures    (from   Ref.     11)    19 

2.5  Effect   of    Saturation    Temperature    on    R-114-0il 
Mixtures    (from    Ref.    6)       22 

2.6  Reentrant   Cavity  Geometry  Factors    ........    k    25 

2.7  Surface    Details  of   Gewa-T   and    Thermoexcel-E 
Reentrant   Surfaces      26 

2.8  Comparison    of    R-11    Boiling    Performance  from 

Various   Commercial   Surfaces     (from    Ref.    4)     28 

Schematic    of    Experimental   Apparatus    31 

Photograph    of    Experimental    Apparatus      32 

Photographs    of    Computer-Controlled-Valve 

Equipment 33 

Oil   Sampling    Apparatus 36 

Sectional    Views  of   Short    Boiling   Tube    37 

Sectional    View    cf    Long    Boiling    Tube 33 

Long    Tube    Thermocouple   Channels 41 

Photographs    of    Boiling   and    Foaming    at    0   and    1 

Percent  Oil 50 

4.2         Photographs    of    Eoiling    and    Foaming    at    3   Percent 

Oil 51 


3. 

1 

3. 

2 

3. 

3 

3. 

4 

3. 

5 

3. 

6 

3. 

7 

4. 

1 

4. 3        Photographs   of    Boiling    and   Foaming    at    10 

Percent   Oil - 52 

5.1  Time   Dependency   of   the   Measured   Heat-Transfer 
Coefficient    in    Press-Fit    Tubes      56 

5.2  Comparison   of    Solid  and    Soft-Soldered   Short 

Tube   Heat-Transfer   Coefficients   58 

5.3  Axial  Temperature    Variation   of    Smooth   Tube      ....    59 

5.4  Axial  Temperature   Variation   of    High    Flux   Tube    ...    60 

5.5  System  Response   to   a   Transient   with 
Computer-Contrclled   Valve   63 

6.1  Heat-Transfer    Performance  of   Smooth   Tube      .....    65 

6.2  Eoiling    Heat-Transfer   Coefficient    of    Smooth 

Tube *. 6  7 

6.3  Relative   Effect  of   Oil   on   Smooth    Tube    Boiling 
Heat-Transfer   Performance    69 

6.4  Heat-Transfer    Performance  of   High    Flux   Surface      .    .    70 

6.5  Boiling    Heat-Transfer    Coefficient    of    High    Flux 
Surface ...72 

6.6  Relative    Effect   of   Oil   on   High   Flux    Boiling 
Heat-Transfer    Performance    74 

6.7  Comparison    of    High   Flux    and   Smooth   Tube  Boiling 
Heat-Transfer   Performance    75 

6.8  Relative    Improvement   of    High    Flax      Surface    over 
Smooth  Tube  versus    Oil 77 

6.9  Effect  of    Saturation   Temperature    on    Boiling 
Heat-Transfer    Coefficient    78 

6.10  Effect   of   Oil    en   Variation   of    3oiling 
Heat-Transfer    Coefficient   with    Saturation 
Temperature 79 

A. 1         Thermocouple   Discrepancy    Correction    86 


I.     INTRODUCTION 

A.       BACKGEOOHD 

The  U.S.  Navy  currently  uses  refrigerant  R— 114  in 
centrifugal  chilled-water  air-conditioning  plants  aboard 
submarines  and  surface  ships.  The  Navy  hopes  to  reduce  the 
size  of  these  units  and  increase  their  performance  by  using 
enhanced  evaporator  and  condenser  surfaces.  An  experiment 
by  Arai  et  al.  [fief.  1]  produced  a  prototype  200- Ton  E-12 
centrifugal  water,  chiller  that  was  28  percent  shorter  in 
length  and  had  a  50-70  percent  improvement  in  the  overall 
heat-transfer  coefficient  by  employing  the  enhanced  surface 
"Thermoexcel  E"  made  by  the  Hitachi  Company.  Comparisons  of 
various  enhanced  commercial  tubes  by  Yilmaz  and  Westwater 
[Ref.  2],  Marto  and  Lepere  [Eef.  3  ],  and  Carnavos  [fief.  4] 
for  various  refrigerants  other  than  R-114  indicated  that  a 
porous-coating-enhanced  surface,  such  as  Union  Carbide's 
"High  Flux,"  will  exhibit  the  best  boiling  heat- transfer 
performance  in   a    pure   refrigerant. 

The  High  Flux  surface  (see  Figure  1.1)  consists  of  a 
sintered  metallic  matrix  bonded  to  a  metallic  substrate. 
The  surface  is  produced  by  coating  a  smooth  tube  with  a 
binder-solvent  mixture  and  then  applying  a  mixture  of  metal 
and  braze  alloy  powder;  the  tube  is  placed  in  a  furnace  to 
evaporate  the  solvent,  binder,  and  melt  braze  alloy  thus 
forming  a  porous  structure  having  multiple  reentrant  cavi- 
ties   to   enhance    nucleation.       [Eef.    5] 

Since  the  High  Flux  surface  will  be  employed  in  a 
refrigeration  unit  using  an  oil-lubricated,  hermetically- 
sealed,  compressor,  some  amount  of  oil  is  always  present  in 
the    evaporator.         Studies   by    Henrici      and    Hesse    [Bef.    6]   for 
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Schematic    and    Scanning    Electron    Micrograph 
(500x)    of    High    Flux    Surface. 
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smooth    tubes      and   Stephan      [ Ref .    7]    for      the   Gewa-T      surface 
(manufactured    by    the    Wieland    Company)       in    R-114-oil    mixtures 
indicate      that      the    heat-transfer      coefficient      of      enhanced 
surfaces   can  be    significantly   altered   by    oil. 

Experimental  data  showing  the  effect  of  oil  concentra- 
tion on  the  heat-transfer  coefficient  of  the  High  Flux 
surface  in  R-114  are  lacking,  thus  motivating  the  present 
investigation.  This  investigation  was  funded  by  the  Eavid 
W.  Taylor  Naval  Ship  Research  and  Development  Center. 
Details  of  the  experimental  apparatus  are  described  by 
Karasabun  [Ref.  8].  The  smooth  copper  tubes  were  supplied 
by  the  TCieland  Company.  The  High  Flux  coated  copper-nickel 
tubes   were   supplied    by  the   Onion   Carbide   Corporation. 

B.       THESIS    OBJECTIVES 

The    objectives   of   this    thesis   are: 

1.  Take  boiling  data  on  a  smooth  tube  in  R-114  with  and 
without  oil  for  comparison  with  the  data  of  ether 
researchers,  and  to  provide  baseline  data  for  evalu- 
ating  the   boiling   performance   of  the   High    Flux    tube. 

2.  Take  boiling  data  on  a  High  Flux  tube  for  various  oil 
concentrations    (0    to    10    percent    by    mass). 

3.  Study  the  effect  of  saturation  temperature  on  the 
R-114   boiling    behavior. 

4.  Attempt  to  sample  oil  locally  in  the  near  vicinity  of 
a  tube  to  investigate  the  possibility  of  an  oil 
concentration  gradient  around  the  tube  during 
operation. 
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II.    HE VIEW    OF    REFRIGERANT-OIL    MIXTORE    BEHAVIOR 

A.       HOCIEATE    BOILIHG    OF    REFRIGERANT-OIL    MIXTORES    FROM    SMOOTH 
TOBES 

In  1963,  Stephan  [Ref-  9]  published  a  milestone  paper  on 
the  influence  of  oil  on  the  boiling  heat  transfer  of  R-12. 
The  effects  he  noted  have  been  observed  in  most  refrige- 
rants, including  R-114.  In  1972,  Henrici  and  Hesse  [Ref.  6] 
updated  Stephan' s  work  for  R-114-oil  mixtures  boiling  from  a 
smooth  copper  tube.  Figures  2-1  and  2.2  summarize  Henrici 
and    Hesse's  results.  Figure    2.1    shows   that      oil    generally 

lowers  the  heat-transfer  coefficient,  and  that  at  high  heat 
fluxes  and  high  oil  concentrations  (10  percent),  the  effect 
grows  more  pronounced  (slope  decreases).  Figure  2.2  shows 
that  at  some  oil/heat  flux  combinations,  the  heat-transfer 
coefficient  may  actually  be  improved  by  the  addition  of  oil. 
Chongrungreong  and  Sauer  [Bef.  10]  suggest  that  the  heat- 
transfer  behavior  of  refrigerant-oil  mixtures  can  be  attrib- 
uted to  5  major  factors:  1)  the  physical  properties  of  the 
refrigerant-oil  mixture,  2)  the  saturation  temperature  (or 
boiling  pressure) ,  3)  the  tube  diameter,  4)  the  surface 
condition  of  the  tube  (roughness) ,  and  5)  the  hydrostatic 
liquid    head   above    the  tube. 

1  •      £llX§ic  al   Properties 

Refrigerant- oil  mixtures  have  significantly 

different  physical  properties  than  pure  refrigerants. 
Jensen  and  Jackman  [Bef-  11]  report  that  density  and 
specific  heat  behave  ideally  in  refrigerant-oil  mixtures, 
but      that   viscosity      and      surface      tension   do     not.  Ideal 

behavior      of   refrigerant-oil      mixture      density   and      specific 
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heat      does      not    mean    linear       behavior.         The   governing    equa- 
tions   are: 


Po£ 


1    -   C 


(2.1) 


and 


c        -    (1    -  C)    c      „+f!c 
P<*  pre  Po« 


(2.2) 


where  subscripts: 

P=density  1   =  liquid 

C   =   oil    concentration  m    =    mixture 

Cp   =    specific    heat  o    =   oil 

r   =    refrigerant 

Jensen  and  Jackman  report  that  current  ref r igerant- cil 
mixture  viscosity  equations  substantially  underpredicted 
their    experimental      data.  No    predictive    equation      has    been 

suggested  for  the  surface  tension  of  refrigerant-oil 
mixtures,  though  Jensen  and  Jackman  developed  a  correlation 
for    R-113-oil    mixtures. 

Henrici  and  Hesse  [Ref.  6]  experimentally  determined 
the  surface  tension  for  the  R-114-oil  mixtures  that  they 
used    in    their    1971    experiment.  As    shown    in    Fig ure    2.3,    the 

surface  tension  of  the  mixture  first  decreased  up  to  an  oil 
concentration  of  2.5  percent,  and  then  increased  continu- 
ously with  increasing  oil  concentration.  This  type  of  non- 
linear behavior  makes  explaining  the  change  in  heat- tra nsfer 
coefficient  of  refrigerant-oil  mixtures,  due  to  the  changing 
physical    properties    of   these    mixtures,       both    difficult,       and 
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possibly  non-general.  The  behavior  of  refrigerant- oil 
mixtures  is  specific  to  the  particular  mixture  components, 
and  may  be  dependent  upon  the  kind  of  oil  being  used.  Some 
qualitative  consequences  of  adding  oil  to  a  refrigerant, 
however,    can  be    noted. 

The  most  observable  result  of  adding  oil  to  refrige- 
rants is  foaming.  Oil  concentrations  above  1  percent  result 
in  significant  amounts  of  foaming  from  nucleate  boiling. 
The  foam  bubbles  form  because  the  R-114  in  the  R-114-oil 
mixture  is  more  volatile  than  oil  and  vaporizes  first, 
creating  a  gas  bubble  surrounded  by  an  oil-rich  layer  (see 
Figure    2.4).  Since  the   bubbles      are    coated    with      oil   film 

with  a  higher  surface  tension  than  the  bulk  liquid  and  a 
have  lower  density,  they  rise  to  the  top  of  the  liquid. 
Because  of  their  surface  tension,  the  bubbles  build  up  on 
the    liquid    surface    to  produce    a   foam    layer. 

This  foaming  action,  which  is  most  pronounced 
between  1-10  percent  oil  concentration  [Ref.  6],  may  affect 
the  heat  transfer  of  tube  bundles  significantly.  For  single 
tubes,  it  is  the  oil  concentration  gradient  which  would  seem 
to  play  the  major  role,  since  the  foam  rises  away  from  the 
tube  surface  and  could  only  interact  with  the  tube  as  it 
sweeps    by    it   from    the   bottom    to    the    top   of    the    tube. 

The  general  decrease  in  the  heat-transfer  coeffi- 
cient upon  adding  oil  to  pure  refrigerants  (recall  Figure 
2.1)  is  subject  to  many  explanations.  Thome  [Ref-  1**]#  in 
an  extensive  review  of  the  literature,  reports  that  the 
first  explanation  for  the  decrease  in  the  heat- transfer 
coefficient  of  mixtures  was  presented  by  Van  Fijk  et  al .  in 
1956.  The  effect  was  explained  as  being  the  result  of  the 
evaporation  of  the  more- volatile  components,  leaving  an 
oil-rich  layer  with  a  higher  local  boiling  point,  which 
increases  the  amount  of  superheat  required  to  continue 
vaporization        and      bubble        growth,  thus      reducing         the 
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heat-transfer  coefficient.  Stephan  and  Preusser  [fief.  12] 
demonstrated  conclusively  that  the  work  of  formation  of 
bubbles  in  a  mixture  is  greater  than  in  an  equivalent  pure 
fluid-  They  concluded  that  the  mixture  heat-transfer  coef- 
ficient is  lower  than  for  the  equivalent  pure  fluid  because 
of  the  resulting  decrease  in  the  bubble  population.  Any 
variation  from  a  decreasing  mixture  heat-transfer  coefficent 
(like  Henrici  and  Hesse  show  in  Figure  2.2)  is  attributed  by 
Stephan  [fief.  13]  to  the  non-linear  variation  of  the  phys- 
ical properties.  Stephan  proposes  that  the  plot  of  surface 
tension  (Figure  2.3),  along  with  thermal-property  varia- 
tions, accounts  for  the  anomalous  rise  in  the  heat-transfer 
coefficient  between  3-6  percent  oil  concentration. 
Chongrongeong  and  Sauer  [Ref.  10]  state  that  it  is  the  rate 
of  heat  diffusion,  governed  by  the  thermal  properties  of  the 
oil-rich  layer,  that  limits  the  bubble  growth  and  that  the 
surface-tension    effects    are   neglible. 

Thome  [fief.  14]  proposes  that  all  of  the  above 
factors,  as  well  as  the  viscosity  variation,  are  important 
in  explaining  the  rise  in  the  heat-transfer  coefficient  for 
some    refrigerant-oil    mixture    and   heat-flux    combinations. 

In  summary,  all  researchers  agree  that  the  physical 
and  thermal  properties  of  a  refrigerant-oil  mixture  are 
important  factors  in  explaining  the  heat-transfer  behavior 
of    mixtures. 

2-      Saturation    Temperature 

It  has  long  been  noted  that  increased  saturation 
temperature  (i.e.,  increased  boiling  pressure)  increases  the 
boiling        heat-transfer  coefficient        of  surfaces        in 

refrigerant-oil  mixtures-  In  1963,  Stephan  [Eef„  9]  found 
that  at  high  oil  concentrations,  the  heat-transfer  coeffi- 
cient of  refrigerant-oil  mixtures  becomes  constant  with 
respect    to      the    saturation      temperature.         Stephan      proposed 
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that  this  is  because  the  addition  of  oil  to  a  refrigerant 
introduces  a  large  diffusion  resistance,  and  that  since  the 
velocity  of  diffusion  is  almost  independent  of  temperature, 
so  should  the  heat- transfer  coefficient  become  independent 
of  temperature  at  high  oil  concentration.  Figure  2.5  shows 
Henrici  and  Hesse's  data  on  the  effect  of  oil  and  saturation 
temperature  in  R-114.  With  no  oil,  the  effect  of  raising 
the  saturation  temperature  is  seen  to  be  a  rise  in  the  heat- 
transfer  coefficient.  With  oil,  raising  the  saturation 
temperature  is  seen  to  cause  a  slight  drop  in  the  heat- 
transfer  coefficient.  This  effect  has  not  been  explained 
yet. 

3 •   Tube  Diameter 

i 

Cornwell,  Schuller,  and  Einarsson  [Ref.  15]  found 
that  for  smooth  tube  diameters  from  6  mm  to  30  mm,  the 
nucleate  pool  boiling  heat-transfer  coefficient  in  pure 
refrigerants  falls  with  increasing  diameter.  The  effect  of 
tube    dianeter    was  correlated    by: 


Nu    =    C    Re2/3  (2*3) 


where 


Re 


=    q  D 


hfR    M 


M  h      D 

Nu    =  


C   =    150    for    refrigerants 
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Since  equation.  (2.-3)  depends  on  physical  parameters,  the 
effect  of  oil  could  he  significant.  No  known  measurements 
of  the  effect  of  tube  diameter  in  refrigerant-oil  mixtures 
have,    as  yet,    been    performed. 

**-      Surface   Condition 

Several  researchers  [Eef.  16r  17,  18]  have  investi- 
gated the  effects  of  surface  roughness  on  the  heat  transfer 
of  pure  refrigerants.  As  surface  roughness  increases,  the 
heat-transfer  coefficient  was  found  to  increase  due  to 
increased      nucleation.  Nishikawa    [Ref.    18]      reported      the 

effect  of  a  variety  of  surface  roughnesses  in  pure  E-114 
over  a  range  of  pressures  from  0.294  MP a  (42  psi)  to 
2.94    MPa    (420    psi)  . 

5  •      il4£2s  ta  tic    Effect 

The      liquid    column      above      the      boiling    surface      may 
generate      large    static      pressures      which      will   increase      the 
boiling   point.         For    E-114   at    0    °C    (32    <>F)  ,         a    0.3    m    (1    ft) 
liquid   head  will   raise   the   saturation   temperature    about   5   °C 
(9    °F) .  In    large      machines,       this      may      be   a      significant 

effect.  However,  for  small  experimental  apparatuses,  the 
effect   is   negligible. 

B.       NUCLEATE      BOILIHG        OF      REFRIGERANT-OIL         MIXTDRES      FROM 
EHHAHCED    SURFACES 

Webb  £Bef-  19]#  in  an  extensive  review  of  the  evolution 
of  enhanced  surface  geometries,  notes  that  the  ability  of 
roughness  to  improve  nucleate  boiling  performance  has  been 
known  for  over  50  years.  However,  it  was  not  until  1968 
that  the  first  commercial  enchanced  surface  was  patented 
[Ref.  20].  Since  then,  the  number  of  commercial  enhanced 
surfaces   has   dramatically    increased      as   the    understanding    of 
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their  design  and  operation  has  grown.  Of  the  many  possible 
methods      for      heat-transfer      enhancement,  two      areas      are 

currently  being  commercially  developed:  1}  fins  and  surfaces 
with    reentrant    cavities,    and    2)     porous  coatings. 

1  *      Fins  and    Surfaces    with    Reentrant    Cavities 

The  ability  of  surface  abrasion,  open  grooves,  and 
fins  to  improve  the  heat-transfer  performance  of  a  smooth 
surface    was   first      studied    in    the    1930's.  Ths    main    diffi- 

culty with  using  surface  abrasion  to  improve  the  heat- 
transfer  performance  is  that  fouling  of  the  surface 
eventually  returns  the  performance  to  that  of  a  non-abraded 
surface.  Studies  in  the  1950*s  and  1960's  centered  on  fins. 
Recent  comparative  studies  fcr  refrigerants  by  Carnavos 
[Ref.  4]  and  Yilmaz  and  Hestwater  [ Ref .  2]  found  that  fins 
and  grooves  result  in  a  50-100  percent  permament  improvement 
in  the  heat-transfer  performance  compared  to  a  smooth  plain 
tube  in  the  same  refrigerant.  Webb  [Ref.  19],  in  his  liter- 
ature review,  describes  how  researchers  in  the  early  1960*3 
found  methods  to  improve  the  performance  of  fins  by  creating 
reentrant  cavities  on  their  surfaces.  Reentrant  cavities, 
such  as  shown  in  Figure  2.6,  act  as  very  stable  nucleation 
sites  and  thereby  enhance  the  heat-transfer  performance. 
For  a  cavity  to  function  as  a  nucleation  site  and  remain 
active,  even  after  the  surface  is  subcooled,  the  mouth  diam- 
eter (D)  must  fall  within  a  critical  range.  Also,  the 
cavity  must  have  a  reentrant  shape  with  a  maximum  reentrant 
angle  (6).  The  optimum  mouth  diameter  (D)  and  reentrant 
angle    (©)    are    functions    of    the    fluid    properties. 

The  Gewa-T  surface,  patented  in  1979  (manufactured 
by  the  TCieland  Company) ,  and  the  Thermoexcel-E  surface, 
patented  in  1980  (manufactured  by  the  Hitachi  Company) ,  are 
two  commercial  surfaces  which  use  modified  fin  shapes  to 
form    the   necessary    reentrant   cavities.      Figure    2.7    shows   the 
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Figure  2.6    Beentrant  Cavity  Geometry  Factors, 
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(a)   Schematic  cross  section  of  the  Gewa-T  surfac 


Pore 


(I))   Schematic  view  of  the  Thermoexcel-E  surfac 


Figure  2.7    Surface  Details  of  Gewa-T  and 
Theraoexcel-E  Reentrant  Surfaces. 
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details  of  these  tube  surfaces.  Carnavos  [ Ref .  4]  found  that 
in  pure  R-11,  the  Gewa-T  surface  outperformed  a  plain  tube 
by  100-200  percent  and  the  Thermoexcel-E  surface  outper- 
formed a  plain  tube  by  300-400  percent.  Curves  1-5  of 
Figure  2.8  show  the  relative  improvement  in  the  boiling 
heat-transfer  performance  of  R-11  that  can  be  achieved  by 
using  mechanically  produced  reentrant  cavities.  Work 
continues  to  optimize  these  types  of  surfaces  for  the 
various  refrigerants  in  use  commercially.  Both  tubes  have 
been  tested  in  ref rigerant-oil  mixtures  and  did  not  show  a 
significant  decrease  in  performance  [Ref.  1  and  7].  The 
cost  of  these  surfaces  is  not  significantly  higher  than  for 
smooth  plain  tubes,  and  the  performance  improvement  is 
dramatic. 

2  •      Porous   Coatings 

The  second  major  type  of  enhanced  surface  is  the 
porous  boiling  surface.  Webb  [Ref.  19]  details  the  various 
production  improvements  and  coating  variations  that  have 
been  made  to  the  original  1968  patent  by  Milton  of  Union 
Carbide.  The  key  to  the  performance  of  the  porous  coatings 
is  their  small  reentrant  cavities,  which  are  interconnected 
by  substrate  tunnels.  The  particles  used  to  make  the  coat- 
ings are  usually  copper  or  aluminum.  According  to  Vebb , 
researchers  have  found  that  the  critical  variable  is  the 
pore   size  rather      than  the   particle   size.  Large    pores   are 

required  for  fluids  with  high  surface  tension  and  high 
thermal   conductivity.  Small    pores   are   optimum      for   fluids 

with  low  surface  tension  and  low  thermal  conductivity  (like 
refrigerants).  Curve  6  of  Figure  2.8  shows  the  relative 
performance  of  the  High  Flux  surface  to  finned  tubes  and 
mechanically  produced  reentrant  surfaces.  Carnavos  [Ref.  4] 
found  the  High  Flux  surface  to  be  700-800  percent  better 
than    a      smooth    tube      in    R-11.         No      known    studies      have   been 
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published  on  the  performance  of  porous  coatings  in 
refrigerant-oil  mixtures.  Some  studies  of  the  heat-transfer 
performance  of  porous  coatings  in  pure  R-114  and 
refrigerant- oil  mixtures  have  been  made  by  anion  Carbide, 
but    their  results   are  not   found   in   the   open   literature. 
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III.    DESCRIPTION    OF    EXPERIMENTAL    APPARATUS 

A.       OVERALL    APPABATOS 

An  overall  schematic  of  the  experimental  apparatus  is 
shove  in  Figure  3.1,  and  a  photograph  is  shown  in  Figure 
3.2.  Karasabun  [Ref.  8]  describes  the  design,  construction, 
and   operation   of  the   apparatus   in    detail. 

The  apparatus  consists  of  two  Pyrex-glass  tees.  Liquid 
B- 1 1  U  is  boiled  in  glass  tee  (1)  and  is  condensed  in  glass 
tee  (2) .  Gravity  drains  the  condensate  from  the  condenser 
hack       to      the      boiling      section.  A      water-ethylene-glycol 

mixture  at  -17  °C  (1  °F)  is  pumped  through  the  condenser 
cooling  coil  via  a  computer-controlled  valve  (VC)  to 
condense  the  R-1  14  vapor.  The  sump  (7)  that  supplies  tne 
water-ethylene-glycol  mixture  is  cooled  by  a  1/2-Ton,  R-12 
air-conditioning    plant. 

Valve  VC  controls  the  R-114  liquid  temperature  and  pres- 
sure. Figure  3.3  is  a  photograph  of  valve  VC  and  the 
computer-controlled  motor  that  operated  VC.  Opening  VC 
causes  more  R-1  14  to  condense  and  lowers  the  system  pres- 
sure. Also,  it  returns  more  subcooled  liquid  to  the  boiling 
section  which  lowers  the  bulX  liquid  temperature.  Since 
data  at  many  heat  fluxes  was  desired  for  a  constant  tempera- 
ture, it  can  be  seen  that  changing  the  heat  flux  without 
adjusting  VC  would  change  the  system  pressure  and  tempera- 
ture. A  computer-controlled  valve  was  thought  to  be  the 
best  way  to  rapidly  return  the  system  to  the  desired  satura- 
tion temperature  following  a  heat  flux  change.  Sections 
III.E  and  IV. D  describe  in  more  detail  the  computer - 
contrclled  valve  and  the  operation  of  the  system  with  the 
computer-controlled    valve   in    use. 
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Figure  3.1    Schematic  of  Experimental  Apparatus. 
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Figure  3.2   Photograph  of  Experimental  Apparatus, 
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(a)  Computer-Controlled  Valve. 


(b)  Computer  for  Computer-Controlled  Valve. 


Figure  3.3 


Photographs  of  Conpnter-Controlled-Valve 
Equipment. 
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Oil  was  added  to  the  liquid  R- 114  by  draining  it  from  a 
glass  oil  cylinder  (5)  .  The  oil  cylinder  was  refilled  as 
needed  through  valve  7-2  from  the  oil  reservoir  (4)  . 

Two  configurations  of  boiling  tubes  were  tested.  Short 
boiling  tubes  were  tested  to  determine  the  correct  assembly 
procedure  to  obtain  data  on  the  normal  431.8  mm  (17  in.) 
long  boiling  tubes.  The  short  tubes  were  cheaper  to  make, 
thus  more  debugging  attempts  could  be  made  by  testing  them. 
Section  III.C  describes  the  details  of  the  construction  of 
the  short  and  long  test  tubes. 

B.   Oil  SAHPLIHG  APPARATUS 

Following  all  data  runs,  an  attempt  was  made  to  sample 
the  local  oil  concentration  in  the  vicinity  of  a  boiling 
tube.  Figure  3.4  shows  the  oil  sampling  apparatus.  3y 
opening  valves  S-1,  S-2,  and  S-4,  the  probe  line  could  be 
purged,  trapping  a  sample  inside  a  flexible  silicon  tube 
that  was  30.5  mm  (12  in.)  long  with  a  3.16  mm  (1/8  in.) 
inside  diameter  using  pinch  clamps.  3y  weighing  the  sample 
tube  and  then  boiling  off  the  R— 114  leaving  behind  the  oil, 
the  mass  percent  of  oil  in  the  R-114-oil  mixture  was  deter- 
mined. The  mass  fraction  of  oil  was  calculated  by: 


,,     „    „ .      ml  -  ml 

Mass  Fraction  =  (3.1) 

m2  -  ml 


wher  e 


ml  =  mass  of  sample  line 

m2  =  mass  of  sample  line  +  R— 114  ♦■  oil 

m3  =  mass  of  sample  line  +  oil  (after  boiling  off  R-114) 
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A  Precisa  Model    80    electronic    mass  balance    was    used   to    weigh 
the      samples.  The        Precisa      Model      80      is        accurate      to 

±    0.0001    g. 

C.       BCI1IHG    TUBE    CONSTBUCTION 

1.      Short    Tubes 

Figure  3.5  shows  the  design  of  the  short  tubes.  The 
short  tubes  were  15.9  mm  {5/8  in.)  in  outer  diameter, 
12.7  mm  (1/2  in.)  in  inside  diameter,  and  203.2  mm  (8  in.) 
in  length.  The  short  tubes  extended  152.4  mm  (6  in.)  into 
the  liquid  E-114  from  the  left  end  flange.  A  25.4  mm  (1  in.) 
long  epoxy  plug  insulated  the  right  end  of  the  tube.  A  1  mm 
{0.04  in.)  thick  copper  disk  behind  the  epoxy  plug  was  soft 
soldered  in  place  to  act  as  a  pressure  barrier.  The  short 
tubes  were  heated  by  a  500-Eatt  240-Volt  stainless-steel 
cartridge  heater.  The  heater  was  6.35  mm  (1/4  in.)  in  cuter 
diameter   and    101.6    mm    (4   in.)    in   length. 

The  first  short  tube  was  made  from  thick-walled 
copper  tubing.  This  tube  was  solid  oxygen- free,  high 
conductivity  (OFHC)  copper.  Four  1.2  mm  (3/64  in.)  diameter 
holes  were  drilled  into  the  wall  of  this  tube  at  a  diameter 
of  12.7  mm  (1/2  in.)  for  thermocouple  channels.  Since  this 
tube  was  solid,  and  had  no  sleeve  interface,  it  did  not  have 
an  interface  resistance.  Consequently,  it  was  the  reference 
tube  against  which  all  other  tubes  were  compared  to  deter- 
mine   the  amount   of   contact    resistance   they    had. 

Six  other  short  tubes  were  made.  Five  were  made  of 
soft  copper  tubing  and  had  sleeves  inserted  into  the  tube  as 
indicated  in  Table  1.  Soft  soldering  of  the  sleeves  to  the 
tubes  was  determined  to  yield  negligible  contact  resistance 
by  comparison  with  the  solid  tube.  The  last  short  tube  (7) 
was  made  of  90-10  ccpper-nickel  and  was  coated  with  High 
Flux    over   the   active    101.6    mm       (4    in.)       long    section.         This 
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(a)  Sectional   view  of  short  tube. 

A-A 


|  50.3      50.3 


(b)  Thermocouple  sleeve  unwrapped  (at  section  A-A)  to  show 
the  relative  locations  of  the  thermocouple  channels 
(all  dimensions  in  millimeters). 


A-A 


(c)  Left-end  view  of  the  short  tube. 


Figure  3.5    Sectional  Views  of  Short  Boiling  Tube. 
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tube  tested  the  compatibility  of  the  soft-solder-assembly 
method  with  the  copper-nickel  tube  '  and  High  Flux  coating. 
Section  V.A  describes  the  results  of  data  taken  on  the  short 
tubes  and  the  selection  of  the  soft-solder  method  of 
assembly  for   the   long  tubes. 

TABLE    1 
Summary   of  Short    Tube   Assembly    Methods 

Remarks 

Solid,  thick-walled  tube  (reference) 

Slide-fit  (0.005  in.  clearance) 

Slide-fit  (0.002  in.  clearance) 

Press- fit  (0.004  in.  interference! 

Press-fit  (0.006  in.  interference) 
Soft-soldered 
ux   Soft-soldered 


2  .   long.  T  ub e 

figure  3.6  shows  the  design  of  the  long  boiling 
tubes.  These  boiling  tubes  were  15.9  mm  (5/8  in-)  in  outer 
diameter,  12.7  mm  (1/2  in.)  in  inside  diamter  and  431.8  mm 
(17  in.)  in  length.  The  center  203.2  mm  (8  in.)  was  the 
active  test  section.  For  the  copper-nickel  tube,  the  center 
section  was  the  only  portion  of  the  tube  that  was  ccated 
with  High  Flux.  The  remaining  114.3  mm  (4.5  in.)  on  either 
side  of  the  center  section  were  smooth  and  unheated,  and  did 
not  nucleate  under  any  heat  flux  or  oil  condition.  Karasabun 
[Ref.  8]  describes  hew  these  end-surfaces  were  treated  by 
the  data-reduction  program  as  an  extended  fin  from  the 
center   section    and    how   their    heat    loss    was    accounted    for. 

The  center  section  was  heated  by  a  1000-Watt 
240-Volt  stainless-steel  cartidge  heater.  The  heater  was 
6.35  mm  (1/4  in.)  in  outer  diameter  and  203.2  mm  (8  in.)  in 
length.         The    heater    was    surrounded      by    a    copper    sleeve    with 
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eight    1.3   mm      by    1.3    mm    (0.050    in.       by      0.050   in.)  thermo- 

couple channels  in  them.  Figure  3.7  shows  the  details  of 
the  channel  layout.  The  thermocouple  hot  junctions  were 
welded  to  the  sleeve.  Appendix  A  describes  the  calibration 
of  the  thermocouples.  The  channels  were  oriented  to  provide 
both  axial  and  circumferential  readings  of  the  tube  inner 
wall  temperature.  The  sleeve  was  soft  soldered  to  the  tube. 
The  data  on  the  short  soft-soldered  tube  closely  matched  the 
reference  solid- tube  data,  and  the  short  sof t-soldered-tube 
data    matched   similiarly      with    long   tube   data.  The    maximum 

circumferential  wall  temperature  variation  in  the  long 
smooth  tubes  was  0.80  K  (0.31  °F)  at  50  kW/m2  compared  to  a 
solid  tule  circumferential  variation  of  0.34  K  (0.61  °F) . 
Section  V.A  describes  in  more  detail  the  circumferential 
variation  of  temperature  that  resulted  using  the  various 
tube  construction  methods.  Some  axial  temperature  variation 
was  experienced  in  the  long  tubes,  particularly  the  long 
High      Flux      tube.  Non-uniform   heat      generation      from      the 

cartridge  heater  is  believed  to  be  responsible  for  the  axial 
temperature  variation.  Section  V.B  descibes  in  more  detail 
the    long    tube   axial    temperature   variation. 

L.       DATA    ACQUISITION    AND    REDUCTION 

A  Hewlett-Packard  3497A  automatic  data  acguistion/ 
control  unit  was  used  to  read  thermocouple  outputs  and  to 
read  an  analog  signal  representing  the  current  and  voltage 
supplied  to  the  cartridge  heater.  A  Hewlett-Packard  9826A 
computer  unit  was  used  to  control  the  HP-34  97A  and  to 
analyze    and   store    the  data. 

Information  was  entered  through  the  computer  keyboard  to 
prompt  the  HP-3497A  to  automatically  scan  each  channel.  All 
thermocouple  measurements  were  accomplished  by  0.245  mm 
diameter    (30   gage)    copper-constantan     (type-T)    thermocouples. 


40 


Top 


^X- 


4  X   2,6  X 


Bottom 


8x 


"I" 


(a)  View  of  the  boiling  tube  thermocouple  locations  as  seen 
from  the  front  of  the  experimental  apparatus. 
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(b)  Thermocouple  sleeve  unwrapped  (at  section  A-A)  to  show 
the  relative  locations  of  the  thermocouple  channels 
(all  dimensions  in   millimeters). 
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A-A 


RIGHT-END  VIEW 

»  A-A 


3  7 

(c)  End  views  of  the  boiling  tube. 

Figure  3.7    Long  Tube  Thermocouple  Channels 
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A  power  sensing  device,  described  by  Karasabun  [Ref.  8]/ 
converted  the  AC  current  and  voltage  values  supplied  to  the 
cartridge  heater  to  a  0-5  V  analog,  DC  signal  for  scanning 
by  the  HP-3497A.  Table  2  lists  the  channel  allocations  in 
the    HP-3497A. 

Folio-wing  data  acquisition  for  each  point,  results  were 
computed  according  to  the  step-wise  procedure  outlined  by 
Karasabun  [Bef.  8],  and  summarized  in  Appendix  3.  Appendix 
B  also  includes  a  complete  listing  of  the  data-reduction 
program    (DRP2)  . 

TABLE    2 
HP-3  497A   Channel   Allocations 

Channel  Purpose 

25-32  Tube  wall   thermocouples    (location    T(1) 

to    T(8)     in    long    tube,    T(5)-T(3)     not 
used    in   short    tubes) 

33  Liquid    R-114    thermocouple    T  (9) 

34  Backup   liquid    R-114    thermocouple    T(10) 

35  R-114   vapor    thermocouple    T(11) 

36  Sump  thermocouple   TM2) 

62  Cartridge   heater    voltage   analog    signal 

63  Cartridge   heater    current    analog   signal 


E.       CCMPOTER-COHTBOLLED    VALVE 

The  computer-controlled  valve  (VC)  was  a  Whitney, 
screwed  bonnet,  regulating  valve  with  3/8  in.  Swagelok 
fittings.  The  valve  travel  was  10.5  turns  from  full  shut  to 
full  open.  The  valve  handle  was  replaced  by  a  101.6  mm 
(4    in.)       hard      rubber  disk    (seen      in    Figure    3.3)  which    was 

rotated   by   the    motor   pinion    gear. 

The  computer-controlled  motor  was  a  General  Electric 
"Minigear  Motor"  with  a  speed  of  105.7  rpm  and  a  torque  of 
3.39  N-m  (30  lb-in) .  The  motor  direction  was  controlled  by 
two    sets      of   Crydom    solid-state      relays    that  acted      to    open, 
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shut,  .  or  hold  the  motor  depending  on  signals  from  the 
computer.  The  computer  controlled  the  amount  of  time  that 
the  valve  was  moving.  Approximately  3  seconds  were  required 
for   the    motor   to   turn  the   valve  one    turn. 

The  flow  through  the  valve  was  checked  by  a  flow  meter 
and  verified  to  be  approximately  linear  (1  turn  =  10  percent 
flow)  for  the  flows  most  often  used.  The  valve  position  was 
tracked  by  a  1 0- turn  potentiometer,  which  was  connected  to 
the  computer-controlled  motor  by  another  101.6  mm  (4  in.) 
hard  rubber  disk  driven  off  the  motor  pinion  gear.  Hard 
rubber  disks  were  used  instead  of  metal  gears  to  avoid 
damage  to  the  gear  teeth  during  the  program  debugging  stage. 
The  rubber  disks  allowed  sufficient  slipping  when,  for 
example,  the  computer  sent  a  valve-open  signal  even  though 
the   valve  was   fully    open. 

The  E-114  liquid  temperature  input  to  the  computer  was 
provided  by  a  separate  copper-constantan  thermocouple 
installed  in  the  same  liquid  R-114  thermocouple  well  that 
the  HP-3U97A  data  acquistion/control  unit  used.  The  E-114 
ther omocouple  emf  for  the  computer  was  amplified  by  an  Omega 
thermocouple  DC  millivolt  amplifier  before  input  in  to  the 
computer.  The  data  acquisition  system  and  the  computer- 
controlled   valve   system   were   completely    independent   systems. 

The  computer  used  for  the  computer-controlled  valve  was 
an  Octagon  Systems  SYS-2A  microcomputer  with  an  Esprit  I 
terminal  connected  via  an  ES-232C  serial  port.  Appendix  C 
lists    the      control    picgram    used.  The    control      program    was 

written  in  NSC  "Tiny  EASIC.  "  The  control  algorithm  simulates 
a  proportional-integral-derivative  (PID)  controller  to  vary 
the  valve  opening  and  shutting  times.  In  the  final  program, 
limits  were  placed  on  the  numerical  value  of  certain  program 
variables  to  prevent  register  overflow,  jamming  the  valve 
fully  open/shut,  and  to  lessen  the  impact  of  system  noise  on 
the    response   of      the    system.         NSC    Tiny    BASIC      is    limited    to 
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integers  from  -32000  to  32000  and  fractional  numbers  are 
truncated.  The  values  of  the  weighting  factors  (A,B,C)  for 
the  proportional  tern  (E) ,  integral  term  (I) ,  and  derivative 
term     (D)    were   determined    by    trial-and-error. 

The   control   algorithm   simplifies    to    the    following   lines: 

60    Input    recruired    temperature    E 
230    Read   B-1T4    temperature    M 

260    Compute   error    (£) ,    and   change   in   error    (D) 
262    Add    1    each    loop    to  integral   sum    (I)     if    error    is   positive 
264    Add   -1    each    loop    to  integral   sum    (I)    if    error    is   negative 
400    Valve   command    V=  (E/A)  +  (B*D) +  (I/C) 
420    If     ("7>0)    then   open   valve 
430    Otherwise   shut    valve 
540    GCTO    230 

Section  V.C  desrites  in  detail  the  system  response  under 
this  algorithm.  Computer  control  of  valve  VC  was  not  used 
during  the  data  taking  as  originally  planned.  Section  V.C 
discusses   the   reasons   for   this    decision. 


44 


IV.  EXPERIMENTAL  PROCEDURES 

A.   INSTALLATION  OP  TUBE  IN  APPARATUS 

Prior  to  installation  in  the  apparatus,  the  boiling  tube 
surface  was  cleaned  with  a  2  percent  Nital  solution  (to 
remove  surface  oxidation  and  oil) ,  rinsed  with  acetone,  and 
air  dryed.  Short  tubes  tested  with  and  without  the  above 
treatment  showed  no  change  in  the  heat-transfer  perf crmance. 
The  treatment  was  effective  in  removing  the  slight  surface 
oxidation  present  following  soft  soldering  without  changing 
the  smooth  or  High  Flux  tube  performance.  Additionally, 
following  testing  with  oil,  a  treatment  was  needed  to  return 
the  High  Flux  surface  to  the  "no-oil"  condition  for  further 
testing. 

After  installing  the  tube  in  the  glass  tee,  the  appa- 
ratus was  evacuated  to  29  in.  Hg  by  the  portable  mechanical 
vacuum  pump  (6)  shown  in  Figure  3.1.  System  pressure  was 
measured  by  a  Marsh  pressure  gage  (30  in.  Hg  to  150  psi 
range,  ±2.5  in.  Hg  and  ±0.5  psi  accuracy).  The  apparatus 
was  left  at  vacuum  for  2  hours  to  check  for  leaks  prior  to 
each  run.  No  noticable  drop  in  vacuum  was  observed  within 
the  accuracy  of  the  pressure  gage.  Next,  the  system  gage 
pressure  was  raised  to  0.19  MPa  (27  psi),  the  saturation 
pressure  of  R- 1 1 4  at  21  °C  (70  °F) ,  by  opening  valve  V7  (see 
Figure  3.1  for  the  configuration  of  the  valves).  An 
Automatic  Halogen  Leak  Detector,  TIF  5000,  was  used  to  check 
for  E-111  leakage.  The  sensitivity  of  this  detector  is  3 
ppm  minimum  concentration.  After  pressure  equalization  with 
the  F.-114  reservoir  (3),  the  reservoir  drain  valve  (V6)  and 
condenser  return  valve  (V5)  were  opened  to  fill  the  boiling 
tee  with  liquid  R-114. 
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Prior  to  installation,  the  left  end  flange  had  been 
scribe  marked  to  indicate  the  liquid  level  corresponding  to 
1600  cc  (2500  gm)  of  liquid  R-114  at  21  °C  (70  °F)  .  This 
was  the  mass  of  pure  R-114  in  the  apparatus  at  the  beginning 
of  the  data  runs.  The  apparatus  was  now  ready  for  taking 
data . 

B.       GENERAL    OPERATION 

Table  3  lists  the  109  data  runs  accomplished  during  this 
thesis      effort      and      their   purpose.  The      data      runs      were 

numbered  sequentially  and  preceeded  by  a  2-  or  3-letter 
prefix    to   indicate    the    tube    type.       The    tube    prefixes    were: 

SS  =    Short    Solid    Tube 

HF  =    High    Flux    Lona    Tube 

HH  =    wieland    Hard   Copper   Long    Smooth    Tube 

SSF  =   Short    Slide-Fit    Tube 

SPF  =    Short   Press-Fit    Tube 

SST  =    Short    Soft-Soldered   Tube 

SHF  =    Short   High   Flux    Tube 

The  short  tube  runs  consisted  of  6  data  points  at  6 
different  heat  fluxes  (usually  59,  37,  22,  14,  8,  and  5 
kR/m2) .  The    normal     tube      runs   consisted      of      at   least      7 

different  heat  fluxes  with  6  data  readings  at  each  heat  flux 
(usually  98,  61,  37,  22,  14,  8,  and  5  kW/m2)  with  some  addi- 
tional low  heat  fluxes  investigated  for  the  0,  3,  and  10 
percent  oil  cases  to  check  for  the  onset  of  nucleate  boiling 
and    hysteresis    in   the  High   Flux   surface. 

In  all  cases,  the  data  set  was  begun  by  starting  the 
cooling  pump  (8)  and  opening  valve  VC  slightly  to  slowly 
cool  the  liquid  R-114  to  the  desired  saturation  temperature. 
The  R-114  vapor  temperature  was  also  monitored  and  when  both 
liquid  and  vapor  temperatures  were  stable  (usually  after 
30-45  minutes)  ,  the  heat  flux  would  be  established,  the 
saturation  temperature  reestablished,  and  the  data  aguisi- 
tion    unit      was   allowed      to    take      data.       Usually,         the    vapor 
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TABLE  3 
Summary  of  Data  Runs 

Run  Tsat  No. 

No.  (°C)  Pts            Remarks 

KH10  Huns  WH01  to  WH10  taken  by  Karasabun 

SPF11  6.7  3  Debug  new  tube 

SPF12  -2.2  6  Press-fit  tube  4  (0.004  in.  interference) 

SPF13  2.2  6  Effect  of  Tsat 

SPF14  6.7  6  Effect  of  Tsat 

SSF15  -2.2  6  Slide-fit  tube  2  (0.005  in.  clearance) 

SSF16  2.2  6  Effect  of  Tsat 

SSF17  6.7  6  Effect  of  Tsat 

SSF18  -2.2  6  Rotate  tube  90  degrees 

SSF19  -2.2  6  Study  low  heat  flux  error  band 

SSF20  -2.2  6  Study  high  heat  flax  error  band 

SSF21  -2.2  6  Rotate  heater  -90  degrees,  tube  fixed 

SSF22  6.7  6  Shift  heater  and  thermocouples 

SSF23  -2.2  6  Repeatability 

SPF24  -2.2  6  Repeatability 

SPF25  -2.2  6  Repeatability 

SPF26  -2.2  6  Shift  heater  and  thermocouples 

SPF27  -2.2  6  Clean  with  Nital,  acetone,  and  air  dry 

SPF28  6.7  6  Repeatability 

SSF29  -2.2  6  Slide-fit  tube  3  (0.002  in.  clearance) 

SSF30  6.7  6  Effect  of  Tsat 

SPF31  -2.2  6  Repeatability 

SPF32  -2.2  6  Time  Variation  Study  (1  day  later) 

SPF33  2.2  6  Effect  of  Tsat 

SPF34  6.7  6  Effect  of  Tsat 

5535  -2.2  6  Short  Solid  Tube 

5536  -2.2  6  Repeatability 

5537  6.7  6  Effect  of  Tsat 
SPF38  -2.2  6  4  days  later 

SPF39  -2.2  6  Shift  heater  and  thermocouples 

SPF40  -2.2  6  7  days  later 

SPF4  1  -2.2  6  10  days  later 

SS42  -2.2  6  Time  Study  Solid  Tube  (no  effect) 

SST43  -2.2  6  Short  soft- soldered  tube  6 

SST44  -2.2  6  Repeatability 

SHF45  -2.2  4  Debug  short  High  Flux  tube  7 

SHF46  -2.2  5  Study  Tvapor  superheat 

SHF47  -2.2  7  Repeatability 

SHF48  -2.2  8  Repeatability 

SHF49  -2.2  6  Study  hydostatic  head  (+1  in.  level) 

SHF50  -2.2  6  -Increase  data  pts  at  heat  fluxes 

SHF51  -2.2  6  -Runs  49-54  form  data  set 

SHF52  -2.2  6  -Each  run  at  different  heat  flux 

SHF53  -2.2  6  -Compare  with  set  55-60 

SHF54  -2.2  6  End  of  Set  49-54 

SHF55  -2.2  6  Study  hydrostatic  head  (+0.5  in.  level) 

SHF56  -2.2  6  -same  as  above  set  49-54 

SHF57  -2.2  6  -same 

SHF58  -2.2  6  -same 

SHF59  -2.2  6  -same 

SHF60  -2.2  6  End  of  Set  55-60 

SHF6  1  -2.2  6  Install  Tvapor  radiation  shiell 

SHF62  -2.2  6  Study  Tliguid  subcooling 

SHF63  -2.2  6  Clean  High  Flux  w/Nital  and  acetone 

SHF64  -2.2  6  Repeatability 

SHF65  -2.2  6  Repeatability 

SHF66  -2.2  6  Repeatability 

SHF67  -2.2  7  Rotate  tube  +90  and  +180  degrees 

SHF68  -2.2  6  Repeatability 
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Table   3 


Summary  of  Data    Runs     (cont'd) 


Run 
No. 

SHF6  9 

SHF70 

SHF7  1 

SHF72 

SHF73 

SHF74 

SHF75 

SHF76 

SHF77 

TCH78 

WH7  9 

FH80 

WH81 

WH8  2 

tfH83 

J7H84 

WH85 

WH86 

WH87 

WH88 

WH89 

WH90 

WH91 

WH92 

WH93 

WH9  4 

■SH95 

WH96 

T7H97 

>TH9  8 

HH9  9 

HF100 

HF101 

HF10  2 

HF103 

HF10  4 

HF10  5 

HF106 

HF107 

HF10  8 

HF109 

HF110 

HF11  1 

HF112 

HF113 

HF114 

HF115 

HF116 

HF117 

HF118 

SPF1  19 


Tsat 
(°C) 

-2.2 

6.7 

-2.2 

-2.2 

-2.2 

-2.2 

-2.2 

-2.2 

-2.2 

-2.2 

-2.2 

-2.2 

-2.2 

-2.  2 

-2.2 

6.7 

-2.2 

-2.2 

6.7 

-2.2 

6.7 

-2.  2 

-2.2 

6.7 

6.7 

-2.2 

-2.2 

6.7 

-2.2 

-2.2 

6.7 

-2.2 

-2.2 

-2.2 

6.7 

6.7 

6.7 

-2.2 

6.7 

6.7 

-2.2 

6.7 

-2.2 

-2.  2 

-2.2 

6.6 

-2.2 

-2.2 

6.7 

6.7 

-2.2 


No. 
Pts 

7 
36 
36 
36 
36 
36 
36 
36 
36 

9 
42 
1  1 
42 
42 
38 
49 
66 
48 
48 
48 
48 
36 
48 
48 
43 
48 
61 
53 
66 
45 
66 
66 
72 
24 
66 
69 
56 
48 
48 
43 
48 
48 
39 
72 
48 
48 
48 
72 
72 
48 
36 


Remarks 


Apply  thermal 
Effect  of  Tsa 
Repeatability 
Debug  oil  add. 
cl 
oi 
oi 
oi 
oi 


grease  to  heater  (no  effect) 


Remove  and 

2  percent 

3  percent 
6    percent 

1 0  percent 
Axial/Circumf 

0  percent  oi 
Rotated  tube 
Repeatability 

0  percent  oi 
Repeatability 

0  percent  oi 

0  percent  oi 

1  percent  oi 

1  percent  oi 

2  percent  oi 

2  percent  oi 

3  percent  oi 
3  percent  oi 
3  percent  oi 
6  percent  oi 
6  percent  oi 

10  percent  oi 

10  percent  oi 

10  percent  oi 

Repeatability 

10  percent  oi 

0  percent  oi 

0  percent  oi 

Repeatability 

0  percent  oi 

0  percent  oi 
Repeatability 

1  percent  oi 

1  percent  oi 

2  percent  oi 

2  percent  oi 

3  percent  oi 
3  percent  oi 
3  percent  oi 
6  percent  oi 
6  percent  oi 

10  percent  oi 
1 0  percent  oi 
10  percent  oi 
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temperature  would  read  higher,  up  to  2  K,  than  the  liquid 
temperature  due  to  the  vapor  .becoming  superheated.  The 
apparatus,  though  coded  by  the  R-114,  was  still  hot  enough 
(60  °F)  to  superheat  the  vapor.  Measurements  of  the  liquid 
temperature,  system  pressure,  and  vapor  temperature  (with 
vapor  probe  shielded  by  a  radiation  shield)  confirmed  this. 
The  liguid  temperature  best  indicated  the  saturation 
temperature. 

The  HP-3497A  data  aguisition  unit  would  scan  each 
channel,  compute  the  heat- transfer  coefficient,  print  the 
results  (an  example  printout  is  shown  in  Appendix  D)  ,  and 
store  the  data  on  the  floppy  disk. 

Following  the  taking  of  all  data  points,  the  data  set 
was  statistically  analyzed  by  subroutine  STATS  to  compute 
the  average  heat-transfer  coefficient  at  each  heat  flux  and 
the  standard  deviation  of  the  6  data  points  for  a  given  heat 
flux.  The  standard  deviation  of  the  6  data  points  was 
usually  0.5  percent  for  the  heat  flux  and  1-2  percent  for 
the  heat-transfer  coefficient. 

After  all  data  sets  at  a  given  oil  concentration  were 
complete,  oil  was  added  via  valve  VI.  The  oil  immediately 
dissolved  in  the  R-114.  No  carryover  to  the  condenser  was 
noted  except  for  several  small  drops  at  10  percent  oil  and 
the  highest  heat  flux  during  the  last  few  data  sets. 
Foaming  occurred  with  the  addition  of  oil,  and  increased 
with  both  increasing  heat  flux  and  increasing  oil  concentra- 
tion. Figure  4. 1  to  4.3  show  photographs  at  heat  fluxes  of 
30  kW/m2  and  98  kW/m2  for  oil  concentrations  of  0,  3,  and  10 
percent.  When  oil  addition  was  not  taking  place,  the  oil 
cylinder  and  reservoir  were  isolated  from  the  apparatus  by 
valves  V3  and  V4  to  minimize  R-114  absorption  by  the  oil. 
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(a)  0  percent  oil  at  100  kW/m  .   No  Foaming. 
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(b)  1  percent  oil  at  100  kW/m  .   Foaming  begins  to  appear. 

Figure  4. 1    Photographs  of  Boiling  and  Foaming 
at  0  and  1  Percent  Oil. 
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(a)  3  percent  oil  at  37  kW/m2. 


(b)  3  percent  oil  at  100  kW/m2. 

Figure  4.2    Photographs  of  Boiling  and  Foaminq 

at  3  Percent  Oil.  3 


5.1 


(a)  10  percent  oil  at  37  kW/m2 


(b)  10  percent  oil  at  100  kW/m2. 

Figure  4.3    Photographs  of  Boiling  and  Foaainq 

at  10  Percent  Oil. 
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C.       SABPLIHG   OIL 

Oil  sampling  was  done  at  room  temperature  to  prevent  the 
temperature  of  the  sampling  lines  from  boiling  off  the  R-114 
and  yielding  false  readings.  The  sampling  procedure  used 
was : 

1.  Weigh  the  empty  flexible  sample  tubing  and  pinch 
clamps.    Connect   them    to   the   sampling   apparatus. 

2.  Open  S- 1 ,  S-4,  and  set  S-2  to  position  2  to  provide  a 
purge  path   for  the   sample   flow. 

3.  lower  the  sample  container  to  a  height  below  the 
level   of    the   R-114   glass   tee. 

4.  After  purging  the  sample  probe  and  lines,  use  pinch 
clamps  to  isolate  the  304.8  mm  (12  in.)  flexible 
sample  tubing  and  trap  a  sample.  Shut  S-1  and  switch 
S-2    to   position    1. 

5.  Disconnect  the  flexible  sample  tubing  and  immerse  it 
in  ice  to  lower  its  internal  pressure  below  atmos- 
pheric  pressure. 

6.  Open  S-3  to  provide  a  vent  for  the  purge  container 
and  pour  the  purge  volume  back  into  the  boiling  tee 
by  lifting  it  to  a  height  so  that  it  flows  by 
gravity. 

7.  Weigh  the  sample  line  (the  mass  balance  used  was 
accurate    to   ±    0-0001    g)  . 

8.  Open  a  pinch  clamp  and  allow  air  to  warm  the  sample 
line  and  evaporate  off  the  liquid  R-114.  Keep  the 
sample  line  on  the  mass  balance  to  collect  any  drops 
of  oil  that  may  splatter  as  the  liquid  3-114  evapo- 
rates. 

9.  After  allowing  several  hours  for  the  R-114  to  evapo- 
rate,   reweigh    the   flexible   sample    tubing. 

10.  Calculate      the   mass      percent      of      oil    using      equation 

(3.1). 
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7.     SYSTEM    OPEBATIPI    AND    PROBLEMS 

A.       CIRCUMFERENTIAL    WALL    TEMPEEATOEE    VARIATION 

Following  the  construction  of  the  solid,  thick-walled, 
reference  tube,  an  investigation  was  begun  to  develop  an 
assembly  method  which  would  result  in  negligible  contact 
resistance  between  the  inner  tube  wall  and  the  copper 
sleeve,  and  yield  data  comparable  to  that  for  the  solid 
tube. 

Short  slide-fit  tubes  2  and  3  (see  Table  1)  had  sleeves 
that  could  be  slid  into  the  test  tube  with  tube  3  having  a 
tighter  clearance.  The  slide-fit  tubes  exhibited  a  circum- 
ferential wall  temperature  variation  of  7-14  K  at  50  kW/m2 
in  pure  R- 1 1 4  boiling  at  -2.2  °C.  This  variation  matched 
data  on  a  similar  long  tube  (with  slide-fit  sleeve)  tested 
by  Karasabun  [Ref.  8].  Tube  2  had  an  average  wall  tempera- 
ture of  33  °C  and  tube  3  had  an  average  wall  temperature  of 
25  °C  at  50  kW/m2.  The  tighter  clearance  of  tube  3  resulted 
in  a  lower  contact  resistance,  and  the  wall  temperature 
subseguently    dropped-  The   short   solid    tube      under    similar 

conditions  exhibited  a  0.34  K  circumferential  wall  tempera- 
ture variation  and  an  average  wall  temperature  of  only 
10.7  °C .  Sauer  et  al.  [Ref-  21 ],  with  a  similar  experi- 
mental apparatus,  used  a  mechanically-press-fitted  brass 
sleeve  for  obtaining  data.  Tubes  4  and  5  had  copper  sleeves 
that    were      mechanically   cold      pressed    into      the    tube.  The 

diametral  interference  was  0.0  1  mm  (0.0004  in.)  for  tube  4 
and  0.015  mm  (0.0006  in.)  for  tube  5.  The  interface  pres- 
sure obtained  was  calculated  to  be  15.2  ?1Pa  (2200  psi)  for 
tube    4    and   22.8    MPa     (3300    psi)     for    tube    5.       During    pressing, 

the      sleeves      were      lubricated      with      glycerin       (C    H  0    )         a 

J     8     3 
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high-thermal-conduc tivity  compound  that,  according  to 
Incropera  and  Dewitt  [Bef.  22],  should  result  in  a  10-100 
times  reduction  in  the  contact  resistance  when  combined  with 
the    contact   pressures  mentioned   above. 

Tubes  4  and  5  had  circumferential  wall  temperature  vari- 
ations of  0.9-1.5  K  with  an  average  wall  temperature  of 
16.2  °C  for  tube  4  and  14.7  °C  for  tube  5  during  initial 
testing.  However,  with  time,  the  circumferential  wall 
temperature  variation  grew  to  5-8  K  and  the  average  wall 
temperature  increased  to  about  28  °C  for  both  tubes.  The 
resulting  drop  in  heat-transfer  performance  is  shown  in 
Figure  5.1.  This  result  is  believed  to  be  due  to  the  phenom- 
enon   referred   to  as   "stress    relaxation." 

Stress      relaxation   is      a      form      of  creep.  The      Metals 

Handbook  [fief.  23]  notes  that  copper  alloys  easily  undergo  a 
decrease  in  stress  resulting  from  transformation  of  elastic 
strain   into      plastic    strain      in  a      constrained   solid.  The 

phenomena  occurs  even  at  relatively  low  operating  tempera- 
tures (at  25  °C,  an  80  percent  drop  in  stress  can  occur  in 
200  hours) ,  and  is  of  most  concern  in  applications  like 
press-fits  and  solderless  wrapped  copper  connectors.  The 
time-dependent  data  of  Figure  5.1  appears  to  be  the  result 
of  a  dropping  press- fit  interface  pressure  yielding  a  higher 
contact  resistance,  higher  inner  wall  temperatures,  and 
lower   heat-transfer    coefficients. 

Stephan  and  Mitrovic  [Bef-  7]  used  a  combination  of 
mechanically  press-fit  and  soft-soldered  inner  sleeves  to 
obtain  data  on  the  Gewa-T  surface  in  R-114.  Soft-soldering 
was  initially  overlooked  by  this  experimenter  because  the 
normal  heating  method  used  is  an  oxy-acetylene  torch  on  the 
tube    surface.  The   flame   temperatures   of      an   oxy-acetylene 

torch  are  above  800  °C  (1500  °F) .  The  resulting  oxidation 
and  heat  damage  (the  High  Flux  coating  melts  at  800  °C)  to 
the      High      Flux      surface   was      unacceptable.  Tests      showed 
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however,  that  the  inner  cartridge  heater  could  be  used  to 
apply  a  controlled  amount  of  heat  to  melt  a  low  temperature 
(melting  point  160  °C)  solder,  while  closely  monitoring  the 
tube      temperature      using    the      sleeve      thermocouples.  This 

method  was  used  to  produce  soft-soldered  short  tube  6  which 
showed  circumferential  wall  temperature  variations  of  0. 8  K 
and  an  average  wall  temperature  of  11.5  °C  at  50  kW/m2. 
Figure  5.2  compares  the  results  of  the  soft-soldered  and 
solid  tubes.  The  agreement  between  the  tubes  is  excellent. 
The  long  boiling  tubes  matched  the  short  soft-soldered  data 
very    closely. 

Tests  made  by  rotating  the  test  tubes  90  degrees  and  180 
degrees  showed  that  the  slight  circumferential  wall  tempera- 
ture variation  of  these  tubes  was  due  to  the  surface  charac- 
teristics of  the  boiling  tube  rather  than  due  to  the 
thermocouples.  Thermocouples      located      near      more      active 

nucleation  sites  of  the  smooth  and  High  Flux  surfaces  had 
slightly  lower  local  wall  temperatures.  The  0.8  K  circum- 
ferential wall  temperature  variation  that  resulted  from 
soft-soldering  the  sleeves  of  the  tubes  is  much  smaller  than 
the  11  K  variation  reported  by  Sauer  et  al.  [Bef-  21],  and 
is  about  the  same  as  the  1  K  variation  reported  by  Stephan 
and    Mitrovic  [fief.    7]. 

B.       AXIA1    TEMPEE2TUBE   VABIATION   OF    LONG    TUBES 

Karasabun  [Bef.  8]  reported  an  axial  temperature  varia- 
tion of  about  20  K  along  the  inner  wall  of  his  slide-fitted 
sleeve.  Stephan  and  Mitrovic  [Bef.  7]  reported  an  axial 
temperature  variation  of  1  K  along  the  sleeve  they  used. 
The  long  soft-soldered  tubes  tested  in  this  experiment 
exhibited  an  axial  temperature  distribution  that  varied  with 
heat  flux.  Figures  5.3  and  5.4  show  the  axial  temperature 
distribution  as      a    function      of    position      along    the      boiling 
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Axial  Location  (cm) 
Axial  Temperature  Variation  of  Smooth  Tube. 
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surface  for  the  smooth  tube  and  High  Flux  tube  respectively. 
The  variation  at  20  kW/m2  (the  heat  flux  at  which  Stephan 
and  Mitrovic  cite  their  1  K  variation)  is  less  than  1  K  for 
both  tubes.  The  maximum  axial  temperature  variation  is  3  K 
for    the    High   Flux    surface  at    98   kW/m2. 

The  cartridge  heaters  used  were  precision- wound , 
ceramic-core,  magnesium-oxide-insulated,  Incoloy-sheathed 
"WATLOfl  FIBEF.OD  heaters.  The  heaters  had  6.35  mm  (0.25  in.) 
long  lava-rock  plugs  at  either  end  of  the  heater  that 
reduced  the  actual  heating  length  on  each  end.  The  heaters 
were  initially  believed  to  generate  a  uniform  heat  flux  at 
all  power  settings,  but  the  axial  wall  temperature  data 
indicate   the   heat   generation   varied    with   the   power    level. 

Since  the  axial  temperature  distribution  was  fairly 
linear  over  most  practical  heat  fluxes  (less  than  37  kW/m2)  , 
the  arithmetic  average  of  all  8  wall  thermocouples  was  used 
to  compute  the  heat-transfer  coefficient.  This  results  in 
slightly  lower  heat-transfer  coefficients  and  is  a  conserva- 
tive estimate  of  the  performance  of  the  High  Flux  surface. 
Appendix  E  analyzes  the  resulting  error  in  the  heat-transfer 
coefficient  from  using  the  arithmetic  average  for  calcu- 
lating  the   heat- transfer   coefficient. 

C.       CCBPUTER-CONTROIIED    VAL7E 

Figure  5,5  shows  the  system  response  to  closing  valve  VC 
manually  and  the  improved  response  of  the  system  when  using 
the  computer-controlled  valve.  3y  correctly  cycling  the 
valve  between  open  and  shut,  the  system  saturation  tempera- 
ture could  rapidly  be  changed-  Establishing  a  stable  equi- 
librium temperature  after  a  large  valve  movement  was, 
however,  difficult  during  the  initial  testing  of  the  appa- 
ratus. The  extensive  trial-and-error  testing  for  the 
correct    weighting    factors   of    the    proportional,    integral,    and 
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derivative  terms  resulted  in  a  large  data  base  concerning 
the  saturated  system  operating  characteristics  and  a  large 
amount  of  practice  for  the  operator  in  manual  control  of 
valve        VC        following  unsuccessful        computer-ccntrolled 

transients. 

Manual  operation  of  valve  VC  was  done  in  conjunction 
with  the  section  of  the  data- reduction  program  (lines  1890 
to  2485  of  program  DEE2  listed  in  Appendix  B)  that  monitored 
the  saturation  temperature  and  the  rate  of  change  of  satura- 
tion temperature  continuously.  It  was  found  that,  because 
of  the  better  temperature  resolution  of  the  HP-3497A  data 
acquisition  unit,  with  manual  operation  and  a  trained  oper- 
ator, temperature  transients  could  be  maintained  within 
±  0.02  K.  This  was  a  much  tighter  control  band  than  was 
possible  with  the  ccmputer-controlled  valve  because  of  the 
poor  resolution  of  its  temperature  sensing  input  from  the 
Omega  thermocouple  DC  millivolt  amplifier.  Additionally, 
after  much  practice  and  experience,  the  manual  method  was 
found  to  be  as  fast  as,  or  faster  than  the  use  of  computer- 
controlled   valve. 

To  obtain  the  highest  accuracy  data  in  the  shortest 
period  of  time,  computer-control  of  valve  VC  was  abandoned. 
Section  VII. B  includes  recommendations  for  further  improve- 
ments to  the  computer  control  system  to  restore  its 
usefulness. 
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(a)  Natural  system  response  to  shutting  valve  VC  to  cause 
a  desired  increase  in  system  temperature. 
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(b)  Improvement  in  system  response  time  using  computer 
control  of  valve  VC  to  reach  a  desired  temperature 
(note  the  different  time  scale). 


Figure  5.5    Systea  Response  to  a  Transient  with 
Computer-Controlled  Valve. 
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VI.  BESULTS  AND  DISCISSION 

A.   BCILIBG  PEBFOBMANCE  OF  SMOOTH  TOBE 

Figure  6.1  shows  the  nucleate  pool-boiling  performance 
of  the  smooth  copper  tune  in  B-114.  The  behavior  with  no 
oil  represents  typical  nucleate  pool-boiling  performance. 
From  point  A  to  point  B,  a  continuous  increase  in  wall 
superheat  ^wo"Tsat^  "*"s  observed  when  heat  flux  is 
increased.  No  bubbles  were  observed  in  this  region  of  the 
curve  as  this  region  represents  natural  convection.  At 
point  B,  incipient  nucleate  boiling  occurs.  From  point  3  to 
point  C  (or  C  for  10  percent  oil),  a  reduction  in  wall 
superheat  is  observed  while  the  heat  flux  is  continuously 
increased.  This  region  is  known  as  the  mixed-boiling 
region,  where  transition  from  natural  convection  to  nucleate 
pool  boiling  takes  place.  In  this  region,  the  heated 
portion  of  the  tube  began  to  activate  an  increasing  number 
of  nucleation  sites,  while  the  unheated  ends  showed  no 
bubbles.  In  fact,  the  unheated  ends  underwent  only  natural 
convection,  due  to  axial  conduction  of  heat  along  the  tube 
wall,  at  all  heat  fluxes.  The  transition  from  natural 
convection  to  nucleate  boiling  occurred  rapidly  when  there 
was  no  oil  present.  The  surface  would  burst  into  nucleate 
boiling  in  less  than  a  second  after  the  first  nucleation 
site  became  active.  At  point  C  (or  C),  all  the  available 
nucleation  sites  were  apparently  active.  After  point  C,  the 
wall  superheat  again  increases  with  increasing  heat  flux  as 
shown  in  region  C-D.  In  region  C-D,  no  new  nucleation  sites 
were  seen  to  become  active.  Instead,  the  bubble  departure 
rate  increased.  When  the  heat  flux  is  decreased  after  having 
established  complete  nucleate  boiling,   the  curve  follows  a 
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different  path  (point  D  to  point  E) .  The  existence  of 
stable  nucleation  sites,  which  remain  active  over  a  wide 
range  of  heat  fluxes,  results  in  better  heat- transfer 
performance  than  in  natural  convection,  resulting  in  a  lower 
wall    superheat. 

The  effect  of  adding  oil  is,  according  to  Stephan 
[Ref.  9],  to  introduce  a  mass  diffusion  resistance  and  lower 
the  heat-transfer  coefficient.  As  seen  in  Figure  6.1,  in 
region  A-B,  both  the  3  and  10  percent  curves  were  lower  than 
the  0  percent  oil  curve.  The  3  percent  oil  curve  is  lower 
than  the  10  percent  oil  curve  probably  because  of  the  non- 
linear physical  property  variations  of  refrigerant-oil 
mixtures.  The  non-linear  variation  of  surface  tension  (see 
Figure  2.3)  would  not  seem  to  be  responsible  for  this 
anomally.  The  curves  in  region  A-B  support  the  contentions 
of  C hongrungeong  and  Sauer  [ Bef .  10]  and  Thome  [2ef.  14] 
that  the  non-linear  variation  of  physical  properties  of 
refrigerant-oil  mixtures,  other  than  surface  tension, 
explains  the  heat-transfer  behavior  of  refrigerant-oil 
mixtures.  The  effect  of  adding  10  percent  oil  was  to  delay 
the  transition  to  ccnplete  nucleate  boiling  on  the  tube. 
With  10  percent  oil,  the  surface  developed  patches  of 
nucleation  sites  that  spread  slowly  with  increasing  heat 
flux,  until  they  covered  the  entire  surface  (point  Cf). 
Region  E-E  (or  Df-E)  shows  that  oil  increased  the  wail 
superheat  slightly  for  3  percent  oil  and  significantly  for 
10  percent  oil.  Again,  this  agrees  with  the  concept  of  an 
increased    mass    diffusion   resistance    by    the    addition   of    oil. 

Figure  6.2  shows  the  heat-transfer  coefficient  of  the 
smooth  tube  in  R-114-oil  mixtures  as  a  function  of  heat 
flux.  The  curves  show  the  heat-transfer  coefficient  of  the 
smooth  tube  in  region  D-E,  after  complete  nucleate  boiling 
has      been   initiated.  The      effect    of      adding      less    than      6 

percent      oil      is      seen      to   be      small       (about      a       10      percent 
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reduction  in  the  heat- transfer  coefficient).  However,  an 
oil  concentration  of  10  percent  causes  a  more  significant 
drop  in  the  heat-transfer  coefficient  (from  0  to  35  percent 
depending   on    the   heat  flux)  - 

Figure  6.3  shows  more  easily  the  degradation  that  oil 
causes  in  the  boiling  heat- transfer  performance  of  the 
smooth    tube.  This    figure    plots   the      heat-transfer    coeffi- 

cient of  the  smooth  tube  relative  to  the  heat-transfer  coef- 
ficient in  pure  R-1  14  as  a  function  of  oil  concentration. 
The  effect  of  oil  can  be  seen  to  depend  also  on  the  heat 
flux.  Oil  can  be  seen  to  generally  degrade  the  performance 
of  the  smooth  tube,  except  at  heat  fluxes  less  than  5  kw/m2 
and      oil  concentrations      between      2      and    3      percent.  This 

behavior  was  also  seen  by  Henrici  and  Hesse  (see  Figure 
2.2).  Since  this  curve  shows  the  heat-transfer  performance 
in  region  D-E,  with  complete  nucleate  boiling,  tne  non- 
linear variation  of  the  physical  properties  of  refrigerant- 
oil  mixtures,  including  surface  tension,  again  probably 
accounts  for  this  anomalous  behavior.  Since  no  measurements 
of  the  physical  properties  of  the  R-114-oil  mixture  were 
made  during  this  investigation,  any  possible  non-linear 
property  variations  of  the  mixture  used  are  unknown, 
requiring   future    work. 

B.       BOILING    PERFORMANCE   OF    HIGH    FLUX    SDRFACE 

Figure  6.4  shows  the  nucleate  pool-boiling  performance 
of  the  High  Flux  surface  in  R-114-oil  mixtures.  The  small 
magnitude  of  the  wall  superheats  obtained  daring  nucleate 
boiling    should    especially   be   noted.  The    High    Flux    surface 

in  pure  R-1 14  showed  typical  natural-convection  region  (A-3) 
behavior.  The  incipient  point  (B)  occurred  at  much  lower 
heat  fluxes  and  superheats  than  for  the  smooth  tube  as  shown 
earlier.  The    transition      to    nucleate      boiling     (B-C)         was 
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similar  to  the  smooth  tube,  occurring  very  rapidly  (less 
than  1  second)  .  Nucleate  boilng  from  the  High  Flux  surface 
results  in  a  drop  of  the  wall  superheat  by  about  a  factor  of 
10.  This  is  due  to  the  extremely  high  density  of  nucleation 
sites  present  on  the  High  Flux  surface.  From  the  low  heat 
fluxes  (1-10  kW/m2)  at  which  the  transition  to  nucleate 
toiling  occurs,  it  can  be  seen  that  the  High  Flux  coating 
also  assists  in  the  activation  of  stable  nucleation  sites. 

Adding  oil  to  the  High  Flux  surface  appears  to  delay  the 
onset  of  nucleate  boiling  (point  B  wall  superheat  increases 
with  increasing  oil  concentration).  However,  the  transition 
to  nucleate  boiling  still  occurred  very  rapidly,  even  at  1 0 
percent  oil.  The  rapid  transition  to  nucleate  boiling  (less 
than  1  second)  on  the  High  Flux  surface  is  probably  due  to 
the  interconnected  cavities  which  can  assist  in  nucleating 
the  entire  surface  once  a  single  site  becomes  initially 
active.  Oil  is  unlikely  to  inhibit  this  characteristic  of 
the  High  Flux  surface  though  it  apparently  delays  the 
initial  activiation  of  the  first  nucleation  site.  As  seen 
in  region  D-E,  for  heat  fluxes  less  than  37  kW/m2,  the 
effect  of  adding  oil  to  the  High  Flux  surface,  once  it  has 
been  nucleating  fully,  is  to  cause  about  a  30  percent 
increase  in  the  wall  superheat.  At  heat  fluxes  in  excess  of 
37  kW/m2  and  at  10  percent  oil,  the  wall  superheat  increased 
dramatically.  As  seen  in  Figure  6.1,  the  wall  superheat  at 
a  heat  flux  of  98  kW/m2  and  10  percent  oil  is  about  the  same 
for  both  the  smooth  tube  and  High  Flux  surface. 

Figure  6.5  shows  the  heat-transfer  coefficient  of  the 
High  Flux  surface  as  a  function  of  heat  flux  for  various  oil 
concentrations.  Again,  oil  is  seen  to  degrade  the  nucleate 
pool-boiling  heat-transfer  performance. 

Air-conditioning  plants  typically  operate  in  the  heat 
flux  range  of  1 0  to  40  kW/m2  with  less  than  1  percent  oil. 
In    this   region    of   pracical    interest,    the    boiling 


71 


m — i — i — i r 


i    1     I I L 


ID 


4-1 

u 

4-> 

4J 

4-> 

4-1 

c 

c 

c 

C 

c 

c 

0) 

<u 

0) 

<u 

01 

QJ 

o 

o 

o 

u 

u 

O 

u 

u 

tH 

u 

u 

u 

u 

o 

<D 

CD 

<u 

<D 

0) 

0) 

a. 

a. 

a, 

Q. 

a. 

CU 

CN 

o 

i— i 

CN 

en 

o 

o 

r— 1 

1 

II     II     II     II     II 


+  oao<ii>  S 


_       t 


01 

E 
\ 


W 


m 


ID 


m 


Q 


a 
V 
•H 
U 
•H 
IH 
*H 

a> 
o 
u 

M 
0)    • 

w  u 

a  ra 

Rj<M 

U  U 
H  => 

I  </1 
4J 

19  M 

D  s 

fflH 

CM 

Cn 

•H  C7> 
•H.H 
•HB5 
O 

mm 
o 


in 

i 

M 

3 
CT> 
•H 


a 


CM 


w  /  /A  >       H 


72 


heat-transfer  coefficient  of  the  High  Flux  surface  is  about 
10  times  better  than  a  smooth  tube.  Also,  it  experiences 
only  a  20  percent  drop  in  the  boiling  heat-transfer  coeffi- 
cient with  the  addition  of  oil  (1  to  10  percent).  This  is 
seen  clearly  in  Figure  6.6  which  plots  the  heat-transfer 
coefficient,  relative  to  the  heat-transfer  coefficient  for  0 
percent  oil,  as  a  function  of  oil  concentration.  From  1  to 
10  percent  oil  at  a  heat  flux  of  14  kW/m2,  the  heat-transfer 
coefficient  is  about  80  percent  of  the  no-oil  heat-transfer 
coefficient. 

Figure  6.6  also  shows  that  the  oil-caused  degradation  of 
performance  on  the  High  Flux  surface  is  nearly  independent 
of   oil   concentration     at    practical   heat   fluxes.  Only   at    a 

heat  fluxof  98  kW/m2  and  6-10  percent  oil,  does  the  High 
Flux  surface  experience  significant  performance  degradation. 
At  high  heat  fluxes  and  oil  concentrations,  the  effect  of 
oil  may  be  to  "clog"  the  interconnecting  cavities  of  the 
High  Flux  surface  due  to  boiling  off  of  the  R-114.  Clogging 
the  R-114  surface  with  oil  would  prevent  replenishment  of 
the  nucleation  sites  with  R-114  liquid,  preventing  the 
nucleation  process  and  leading  to  higher  superheats.  The 
time-dependent  behavior  of  the  High  Flux  surface  in  high  oil 
concentrations  and  at  high  heat  fluxes  was  not  studied  in 
this   experiment. 

C.       CCHPABISON      OF         BIGH      FLOX         TO       SMOOTH         TUBE      BOILING 
PEBFOBMANCE 

Figure  6.7  shows  the  heat-transfer  performance  of  both 
the  High  Flux  surface  and  the  smooth  tube  as  a  function  of 
heat  flux.  Again,  the  7-10  times  improvement  in  the  heat- 
transfer  coefficient  by  the  High  Flux  surface  is  easily 
seen.  At  extremely  high  heat  flux  and  high  oil  combinations, 
the  High  Flux  surface  performs  only  sligthly  better  than  the 
smooth    tube. 
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Figure  6.7  also  shows  that  the  heat-transfer  curves  of 
the  High  Flux  surface  and  the  smooth  tube  are  not  parallel. 
At  low  heat  fluxes,  the  High  Flux  surface  is  more  effective, 
in  comparison  to  the  smooth  tube,  in  enhancing  nucleation 
than  at  moderate  heat  fluxes.  At  high  heat  fluxes,  the 
performance  of  the  High  Flux  coating  begins  to  converge 
toward  the  smooth-tube  curve  because  the  surface  became 
vapor  blanketed.  This  is  expected  since  both  tubes  should 
perform  about  the  same  when  they  are  completely  vapor 
blanketed. 

Figure  6.8  shows  the  relative  improvement  of  the  High 
Flux  surface  over  the  smooth  tube  as  a  function  of  oil 
concentration.  For  the  heat  fluxes  of  practical  interest  in 
air-conditioning  plants,  the  High  Flux  surface  is  7-10  times 
better   than   the    smooth   tube. 

D.        EFFECT    OF    SATURATION    TEMPERATURE    ON    BOILING    PERFORMANCE 

As  reported  by  Stephan  [Ref.  9],  the  effect  of 
increasing  the  saturation  temperature  for  both  the  smooth 
tube  and  High  Flux  surface  was  increased  heat- transfer 
performance.  Figure  6.9  shows  the  improvement  in  heat- 
transfer  performance  in  pure  R-114  achieved  by  raising  the 
saturation      temperature      from      -2.2    °C       (28    °F)  to      6.7    °C 

(44  °F)  .  At  high  heat  fluxes,  little  improvement  is  seen  in 
the  High  Flux  surface  performance  because  the  surface  is 
nearly    vapor   blanketed   with    bubbles. 

Figure  6-10  shows  the  effect  of  saturation  temperature 
on  R-114-oil  mixtures.  Increased  saturation  temperature  is 
again  seen  to  improve  the  performance  of  the  High  Flux 
surface  as  well  as  the  smooth  tube,  even  with  10  percent 
oil.  This   is      consistent    with      the      no-oil   results,         tut 

contradicts   Henrici    and    Hesse's    data     (see    Figure    2.5). 
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E.   10CAI  OIL  SAHPLIHG  EFFOBTS 

The  results  of  the  attempt  to  sample  oil  locally  near 
the  boiling  tube  were  inconclusive.  The  method  outlined  in 
Section  IV. C  did  not  produce  either  repeatable  or  accurate 
results.  Checks  of  the  sampling  method  were  made  by 
sampling  the  bulk  liquid  with  0  and  10  percent  oil  at  zero 
heat  flux.  The  0  percent  oil  check  yielded  oil  concentra- 
tions from  0  to  2  percent.  The  10  percent  oil  check  yielded 
oil  sample  concentrations  from  5  to  25  percent. 

Further  refinements  to  the  sampling  apparatus  are  needed 
to  permit  an  accurate  local  sample  of  oil  in  the  vicinity  of 
a  boiling  tube.  Section  VII. B  contains  recommendations  to 
improve  the  sampling  apparatus. 
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VII •    CONCLUSIONS    AND    RECOMMENDATIONS 

A.       CCHCIDSIOIS 

1.  In  pure  R-114,  the  pool-boiling  heat-transfer  coeffi- 
cient of  the  High  Flux  surface  is  about  10  times 
larger   than   that   of  a    smooth    tube. 

2.  The  High  Flux  surface  began  nucleate  boiling  at  low 
heat  fluxes,  about  1  kW/m2r  compared  to  a  heat  flux 
of  about    10   kfl/m2    for   the   smooth    tube. 

3.  Oil  delayed  the  onset  of  nucleate  boiling  with  the 
High  Flux   surface. 

4.  Oil  resulted  in  about  a  20  percent  reduction  in  the 
heat-transfer  coefficient  of  the  High  Flux  surface 
for  most  practical  heat  flux  (less  than  37  kW/m2)  and 
oil  combinations    (less    than   6    percent)  . 

5.  At  heat  fluxes  of  98  kW/m2  and  greater  than  6  percent 
oil,  the  performance  of  the  High  Flux  surface 
degraded  by  as  much  as  80  percent.  The  performance 
of  the  High  Flux  surface  was  little  better  than  the 
smooth    tube   at  98   kW/m2    and    10   percent  oil. 

6.  The  boiling  heat-transfer  coefficient  of  the  High 
Flux  surface  was  about  7  times  better  than  that  of  a 
smooth  tube  for  oil  concentrations  from  1  to  10 
percent  over  the  range  of  heat  fluxes  employed  in 
air-conditioning   plants    (10-40   kW/m2) . 

7.  The  boiling  heat-transfer  coefficient  of  both  the 
High  Flux  surface  and  the  smooth  tube  increased  with 
increasing  saturation  temperature.  The  improvement 
decreased  at  high  heat  fluxes  for  the  High  Flux 
surface   due    to  vapor   blanketing    of    the   surface. 


81 


8.  Local  oil  sample  results  were  not  accurate  or  repea- 
table  enough  to  determine  if  a  variation  in  the  local 
oil  concentration  occurs  near  a  boiling  tube. 

B.   BECGHHEHDATIOHS 

1.  The  present  cartridge  heaters  should  be  replaced  with 
more  reliable  heaters  that  will  produce  a  uniform 
heat  flux  axially  at  all  heat  loads. 

2.  A  solid  copper  tube  should  be  coated  with  High  Flux 
and  tested.  Cue  to  the  very  small  superheats  of  the 
High  Flux  surface,  even  a  small  amount  of  contact 
resistance  could  affect  the  data  at  low  heat  fluxes. 

3.  The  physical  properties  of  the  R-114-oil  mixtures 
tested  should  be  measured  to  obtain  information  to 
better  explain  the  reasons  for  the  heat- transfer 
performance  of  both  the  High  Flux  surface  and  the 
smooth  tube.  Particularly,  the  anomalous  rise  in 
heat-transfer  performance  of  smooth  tubes  when  1  to  6 
percent  oil  is  present  in  R-114  should  be  studied. 

4.  A  secondary  heater  should  be  installed  in  the  boiling 
section  to  keep  the  total  heat  input  constant.  The 
heat  input  to  the  secondary  heater  should  be 
increased  when  the  heat  input  to  the  boiling  tube  is 
decreased  and  vice  versa.  This  modification  will 
maintain  a  constant  heat  load  on  the  condenser  and 
eliminate  the  need  to  operate  valve  VC,  except  to  set 
the  saturation  temperature  at  the  beginning  of  a  run. 
Alternately,  the  effort  to  computer  control  valve  VC 
should  continue  by  obtaining  an  accurate  thermistor 
with  a  high  temperature  resolution  for  use  with  the 
SYS-2A  microcomputer-  Improved  temperature  resolu- 
tion would  allow  the  computer-controlled  valve  to 
operate  properly  and  free  the  operator  from  the 
demands  of  manual  control  of  valve  VC, 
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5.  The  performance  of  a  bundle  of  High  Flux-coated  tubes 
in  R-114-oil    mixtures    should   be    studied. 

6.  Data  should  be  obtained  on  the  High  Flux  surface  over 
a  wider  range  of  temperatures  and  with  oils  of 
varying    viscosities. 

7.  The  High  Flux  surface  should  be  tested  for  time- 
dependent  heat-transfer  performance  by  taking  data 
periodically    over  a    long    time   during    boiling. 

8.  A  small  oil  sample  container  with  valves  and  a  vent 
line  should  be  manufactured.  The  flexible  tubing  and 
pinch  clamps  used  in  this  experiment  did  not  properly 
hold  or   vent    the   sample. 
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APPENDIX    4 
THEBMOCOOPLE    CALIBRATION 

Karasabun  £Ref.  8]  describes  the  thermocouple  calibra- 
tion equipment  in  detail-  Two  thermocouples  were  calibrated. 
One  was  made  from  wire  at  the  beginning  of  the  roll,  the 
other  from  the  end  of  the  roll,  following  the  making  of  all 
thermocouples   used   in  the   apparatus    and   tests. 

Essentially,  the  manufacturer-supplied  calibration  equa- 
tion for  the  thermocouple  wire,  a  seventh  order  polynomial, 
was  corrected  slightly  by  adding  to  it  a  second  order  curve 
fit  of  the  variation  of  the  manufacturer's  computed  tempera- 
ture for  a  given  emf  from  a  known  set  of  reference  tempera- 
tures (measured  using  a  Hewlett-Packard  2804A  guartz 
thermometer  with  a  temperature  resolution  of  ±  0.0001  K  and 
accuracy   of   ±    0.03    K) . 

The  manufacturer's  emf  to  temperature  conversion  equa- 
tion   is: 

T    =   aQ    +   a-^E      *    a2E2    ♦   a3E3    *  (A.1) 

a    E4    +    a   E5    +    a    E6    +    a    E7 
4  5  6  7 


where 


T  =   temperature    (°C) 

aQ  =0.1 00860910 

ax  =    25727. 94369 

a2  =    -767345.8295 

a3  =   78025595.81 

a4  =   -9247486589 

ac  =   6.97688E+11 

b 

a6  =   -2. 66192E+13 

a?  =   3.94078E+14 

E  =    thermocouple   reading     (volts) 
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Figure  A.1  shows  the  quartz  thermometer  reading  minus 
the  thermocouple  readings  (discrepancy)  versus  temperature. 
The  two  thermocouples  agreed  to  within  0.05  K  of  each  other 
and  the  manufacturer's  seventh-order  polynomial  needed  about 
a  0.  1  K  increase  to  more  accurately  convert  the  emf's  to  the 
true    temperature.  The   correcting      second-order   polynomial 

was: 


DCP   =    b       +    b    T   +    b?T2  (A. 2) 


where 

DCP  =   discrepancy    (K) 

b      =    8.6268968E-2 

0 

b      =   3.7619902E-3 

1 

b      =  -5.0689259E-5 

2 
T      =   thermocouple   reading 

(from   equation    A.  1)     (°C) 

Thus,  the  temperatures  computed  by  the  data-reduction 
program  (DRP2)  were  emf's  converted  to  temperature  by  equa- 
tion A.1  with  corrections  for  that  temperature  computed  by 
equation  A. 2  added  to  the  temperature  to  get  the  true 
temperature. 

Since  the  data-reduction  program  utilized  differences 
between  thermocouples  in  all  computations,  such  as  wall 
temperature  minus  saturation  temperature,  the  corrections 
above  were  necessary  only  for  the  computation  of  items 
dependent  on  the  absolute  temperature,  like  the  fluid  prop- 
erties. Appendix  E  describes  in  detail  the  uncertainty 
analysis  and  the  effect  of  wall  temperature  variation  on  the 
computation  of    the    heat-transfer    coefficient.         Thermocouple 
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variations  of  the  soft-soldered  tubes  were  most  likely  due 
to  a  slight  amount  of  contact  resistance  and  to  the  surface 
characteristics  of  the  tubes  tested,  since  data  runs 
involving  shifting  the  wall  thermocouples  did  not  affect  the 
data. 
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APPENDIX    ? 
DATA    REDOCTION    PROGRAM 

The      data      reduction      program        below      consists      of      the 
following   sections: 

Main  Program  -  Menu   of   subprogram   options 

Sub    Main  -  Take    data    or   reprocess   data 

Sub    Elot  -  Plot   data   on    log-log   scale 

Sub   Poly  -  Compute   least-squares   curve  fit   of   data 

Sub    Plin  -  Plot  data   on  'linear-linear  scale 

Sub    Stats  -  Compute    average   and   standard    deviation    of    data 

Subprogram    Main   consists   of  the    following  steps: 

1.  Create   data    file  for    data   and    plot   points. 

2.  Select    tube    type. 

3.  Monitor  heat  flux  or  saturation  temperature  to  estab- 
lish  steady-state   conditions. 

4.  Scan  all  channels  listed  in  Table  2  and  save  in  data 
file. 

5.  Convert  raw  eif's  to  temperatures,  current,  and 
voltage. 

6.  Compute  the  heat-transfer  rate  for  the  cartridge 
hater. 

7.  Compute  the  average  wall  temperature,  the  wall  super- 
heat,   and   the    film   temperature. 

8.  Compute  the  physical  properties  of  R-114  using  given 
correlations    at    film    temperature. 

9.  Compute  the  natural-convection  heat-transfer  coeffi- 
cient of  R-114  for  the  non-boiling  ends  of  the  test 
tube. 

10.  Compute    the    heat    loss    from    the   non-boiling   ends. 
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11.  Calculate  the  corrected    heat   flux      from    the    test   tube 
to   the   liquid    B-114. 

12.  Calculate      the  boiling      heat-transfer  coefficient     of 
the   B-114    from  the   test   tube. 

13.  Print   data.    Store   heat  flux   and   wall   superheat    values 
in   plot  file. 

The  following  is  a  listing  of  the  complete  data  reduc- 
tion program  (DEP2)  written  in  Basic  3.0  for  the 
Hewlett-Packard   9826    computer. 

I000I  MIT  NAME:  UPF2 

1 005 1  DATE-'       Octcbe-  19,  1 984 

lUlpi  REVISED:    March  18   I  985 

leusi 

I  020  DEEP 

I025  PRINTER  is  I 

I0JQ  FRINf  USING  "4X  "Select  ootion:"'- 

I035  PRINT  USING  'EX  ""0  luting  data  or    r  e-proce  3  3 1  ng  orevious  data""" 

1040  PRINI  USING  "6X.--I  Clotting  data  on  Leg-Log  """                              f 

1045  FfilNI  USING  "6.X,  "-2  Hotting  data  on  Linear""" 

1050  PRINT  USING  -EX,  •■?■  Male  cross-plot  coeffi  file""" 

less    input  Mu 

I0G0       IF     ldo>0    IIICN    CALL    Main 

1065      It-     Icto'l     MIEN   CALL    r  I  „•  t 

1070       ||"     Ido'2     UK  M   COtl     PI  in 

1075       IT     Idu-3    IHFN   CALL    Ceef 

I0P0  EflO 

1095      St'B    Hjin 

1090  COM    'Cc    d  7  )  ,  lc-)l 

1095  DIM   EmfUZ)    r<l2>   Dla(S)   L>2a<6  >  ,Ula<  G  )    Doa<  (i  >  ,La<S  )  ,Lua<  6  )  ,Kcua<fi  ) 

1100  UAIA    0.I0OUG09I    25727.94369      767f<45 .  8295  ,  780:5595  .  8  I 

1105  t'ATti    -924.748G5H9  6.976Bi«?H  I  ,-2.R6l92£»13    3.  940  7BF»  l  4  / 

1110  l.'IV.n    CI  •  ) 

I  I  l'i  I'lMNini    I  3    701 

I  120  fl  r/IH    709 

I  1 25  PI  El' 

1130  IMrul    "ENIFR   MONIH,    DfUE    ANO    1IME    (HM:n0O«H:HH:SS>'  .Date* 

!  135  nt'lI'M!    709;  "  IP"  iCate* 

I  140  OUM'UI  ?09i  "  10" 

I  145  FNirn  709tDale* 

1150  PRINT 

1155  PRIMI     "  Mentli,    itata    jnit    lime     :"i0«tc* 

IIG0  PRIMI 

1165  PRINT  USING  "  I  eX  ,"•  NO  If  :  Program  rune  :  DRP2-"" 

1170  PEEP 

1175  INPUI  "ENTER  DISK  NUHPER*  On 

1180  PRINT  USING  "ISX,""0lsl  nnnber   -  ""  7/"i(Jn 

I  195  PI  £f 

1190  IMrui   ENIFR  INPUT  rft'f'E  I  0'  '054  A  ,  I  =  f  1 1  E  I "  ,  In 

1195  PI  r  r 

1280  INPUI    "CNILR    THERMOCOUPLE    TYPE    (B'NCW ,  I  «01D )"  ,IcaJ  i 

l?05  IF     In-0    THEN 

12  10  PFlP 

1215  INPUI     "GIl'E    ft    NAME    FOR     IMF    RAW    OAT  A    FILE"    D2_file« 

i:Z0  PRINT    USING    "IRX,""New    file    name:     "  *  ,  I  4A"  l02..f  i  le« 

1225  PFTP 

1230  iult'l    "INPUI    SIZE    Of    l  HE    BUflf.Sizel 

i:35  rilfAIF    L'I'AI    02. fl  let    i\  :H 

1240    assign  of  Hi-:  io  n:   lu-t 
1:151 

12501    DL'MMi    FILE    UNTIL    Ni  un    KflOUN 
i:5S      01.  F  1  let-  "DUMMY " 

i:bq    cue/me  hpai  oi.f iu*  .  i  :ei 

;  :C5       rtSSK  N    PI    1  le  I     10    PI     f  1  let 
i:?0      OUIFHT    £f  1  lel  ilJalet 
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1275   BEEP 

i:sa  input  -gh-e  a  name  for  the  plot  file-  p_fiie* 
i2ss  PEEP 

1290   [«PUT  'INPUT  SIZE  OF  FILE  BDAT"  ,Si;e2 

1295   CREATE  BOAT  PJileS  Si?e2 

I 300  ASSIGN  9P)ot  TO  P_ f  t lei 

1305   BEEF 

1310   INPUT  "ENTER  NUMBER  OF  DEFECTIVE  TCS  ( 0'OEFAULT  >"  ,  Idtc 

1315   IF  Idtc-0  THEN 

1 320  Ldtci-0 

I32S  Ldtc2-0 

1330  PRINT  USIN6  MSX,"-No  defective  TCs  exist"-" 

1335  ENO  IF 

1340  IF  Idtc-I  THEN 

1345  BEEP 

1350  INPUT  "ENTER  DEFECTIVE  TC  LOCATION"  ,Ldtc 1 

1355  PRINT  USING  "I6X,""TC  is  defective  at  location  " "  ,0 "  i  L dt  c I 

I3S0  Ldtc2«0 

I3SS  ENO    IF 

1370  IF  Idtc-Z  THEN 

1375  BEEP 

1330  INPUT  'ENTER  DEFECTIVE  TC  LOCATIONS"  .Ldtc I  ,Ldtc2 

I3BS  PRINT  USING  ■|GX,""TC  are  defective  at  locations  " "  .0  ,4X  ,D"  sLdtc I  ,Ldtc2 

1390  ENO  IF 

1395  IF  Idtc2  THEN 

1400  BEEP 

1405  PRINTER  IS  I 

1410  BEEP 

ui5  rRinr  "invalid  entry- 

1420  FRINTER  IS  701 

1425  GOTO  1305 

1430  ENO  IF 

1435  OUTPUT  SFiIel  iLdtcl  Ldtc2 

I  4401  [n»|  oo I  ton 

1445  ELSE 

1450  BEEP 

1455  INPUT  "GIVE  THE  NAME  OF  THE  EXISTING  OAfA  F ILE "  ,D2_f l 1 eS 

I4G0  PRINT  USING  "IBX,""Old  file  name:  " "  ,  I 4A" 1 02_f 1 leS 

I4G5  ASSIGN  @File2  TO  02_fileS 

I  470  ENTER  9F  1 le2iNrun 

1475  ENTER  @F  i  le2 ;OoldS 

1400  BEEP 

U85  rNPUT  "GIVE  A  NAME  FOR  PLOT  FILE"  ,P_f  i le* 

1490  8FEP 

1495  IMFUT  "INPUT  SIZE  OF  FILE  80AT".Size2 

1500  CREATE  BDAT  P_f  1  lei  Si ;e2 

1505  ASSIGN  SPIot  TO  P..file$ 

1510  PRINT  USING  "IBX.'Thls  data  set  talen  on  1  " '  ,  UA"  iDoldi 

1515  ENTER  ?r  l  le2  iL  dt  c  I  ,Ldt  c2 

1520  IF  Ldtcl  0  OR  Ldtc2'0  THEN 

1525  TRINT  USING  '  I  6X  .  "  "  Therixocouo  les  u/ere  defective  at  loca  t  Ions  :  • "  ,2(  30  ,  4*  ) " 

Ldtcl ,Ldtc2 

1530  ENO  IF 

1535  ENTER  9Fi le2il  tt 

1540  ENO  IF 

1545  PRINTER  I?  I 

1550  IF  ln  =  0  THfN 

1555  BEEP 

1560  PRINT  USING  '4X,  ""Select  tube  type""" 

I5S5  FRINT  USING  "GX  ,  " "3=Snocth  4  inch  Ref""" 

1570  PRINT  USING  "SX  .*" I 'Smooth  4  inch  Cu  (Press/Slide)""" 

1575  PRINT  USING  "6X,""2*Soft  Solder-  4  inch  Cu""" 

1580  PRINT  USING  "6X,"*3-Sof»  Solder  4  inch  HIGH  FLUX""" 

1585  PRINT  USING  "BX , "  "4-Wieland  Hard  8  inch"" 

I5')0  PRINT  USING  "SX  .""S-HICH  FLUX  8  inch""" 

1595  PRINT  USING  "6X  , "  "6-6EUA-K  19  F'm""" 

1600  INPUT  lit 

1605  IF  I  t  t  6  THEN 

1610  BEEP 

1615  PRINT  "INVALID  ENTRY" 

1620  GOTO  I5R0 

I62S  END  IF 

1630  OUTPUT  aFilel i It t 

IB35  ENO  IF 

1610  FRINTER  IS  701 

IG45  PRINT  USING  " I 6X  .  ' " Tube  Tvoa:      ""  0" i 1 t t 
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1650  BEEP 

1655  INPUT  "ENTER  OUTPUT  VERSION  < 0=LONG  ,  I -SHORT  2-NONE > "  , Iov 

16601 

1 665 !  OI*Oianeter  at  thermocouple  oosltions 

1670  DATA  .01 II 125,  .01 II 125  ,  .011  I  I  25  ,. 01 29540  ,. 01 2  446 , .0129540, .0100965 

1675  REAO  Dla(  •  > 

I6P0  DI-Ola( Itt  I 

16851 

1690'  02*Olameter  of  test  section  to  the  base  of  fins 

1 595  OATA  .015875  ,  .0I587S  ,  . 015875  ,. 01 5824  ,  .015875.  .815824  ,  .01270 

1700  REAO  02a( •  ) 

1705  02-02a( Itt  > 

1710! 

1715'  Di-Inside  diameter  of  unenhanced  ends 

1720  DATA  .0127, .0127  ,  .0127  ,  .0132  .  .0127  .  .0132  ,  .01  I  I  125 

1725  READ  Dia( •  ) 

1730  Di=0ia<  Itt  ) 

1735! 

17401  Do-Outside  diameter  of  unenhanced  ends 

1745  OATA  .015875.  .015375. . 01 5875  ,  . 015B24  ,  .015875. . 0  I  5824 ,. 01270 

1750  REAO  0oa( •  ) 

1755  Do=Doa(  Itt  ) 

17601 

1765'  L-Lenoth  of  enhanced  surface 

1770  OATA  . 1016  .. 1016  ,. 1016  .. 1016  ,  .2032  .  .2032  ,  .2832 

1775  REAO  Lai  •  I 

1780  L-LaUtt) 

1785'  , 

17901  Lu-Lenglh  of  unenhanced  surface  at  the  ends 

1795  DATA  .0254 , .0254 , .0354 ,.0254 , ,0762 , .0762 . .07B2 

1800  REAO  Lua( •  ) 

1805  Lu-Lua< Itt  > 

1810' 

1815'  Kcu= Therma L  Conductivity  of  tube 

1820  DATA  398  344,344,45.344,45  344 

1825  REAO  Kcua( • I 

1830  Kcu«Kcua< Itt) 

1835  A-PI«'0o  2-Di -2 >/4 

1840  P-PI»Do 

1845  J-l 

16^0  S--0 

1855  Sy-0 

I860  S<s  =  0 

1865  S.-y  =  0 

1870  Recieat  :  ' 

1875  If"  lm=0  THEN 

1880  ON  KEY  0,15  RECOVER  1870 

1885  FRINTER  IS  I 

1890  PRINT  USING  "  4X  ,""  SELECT  OPTION"" 

1895  PRINT  USING  "6X  , " "0* T AKE  DATA""" 

IS00  PRINT  USING  'BX,"1-SET  MEAT  FLUX""" 

1905  PRINT  USING  "6X,""2=SET  Tsai'"" 

1910  PRINT  USING  " 4X  , " "NOTE :  KEY  0  »  ESCAPE""" 

1915  BEEP 

1920  INPUT  Ido 

1925  IF  ldo=0  THEN  2495 

1930' 

1935'  LOOP  TO  SET  HEAT  FLUX 

1940  IF  Ido-I  THEN 

1945  OUTPUT  709, "AR  AFG2  ALG3  UR5" 

1950  BEEP 

1955  INPUT    "ENTER    OESIRED   QclD-.Dado 

I960  PRINT    USING    "  4X  , "  "DES  IRED    Oflo       ACTUAL    Odo"" 

1965  Err- 1000 

1970  FOR     1=1     TO    2 

I97S  OUTPUT    709i"AS    SA" 
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1 980  Sum-0 

1985  FOR  Jl-I  TO  5 

1990  ENTER  709i£ 

1995  Sum-Sun+E 

2000  NEXT  Ji 

2005  IF  I-l  THEN  Volt-Sum'5 

2010  IF  1-2  THEN  Anp=E 

2015  NEXT  I 

2020  Atjdp-VoH«Aiip/'(PI«02«L> 

2025  IF  ABSfAqdo-Oqdg >  Err  THEN 

2030  IF  AqdD  Dadp  THEN 

2035  BEEP  4000  .2 

2040  SEEP  4000, .2 

2045  BEEP  4000  .2 

2050  ELSE 

2055  BFEP  250  .2 

2060  BEEP  250,. 2 

2065  BEEP  250,. 2 

2070  END  IF 

2075  PRINT  USING  " 4X  ,  MZ . 30E  ,  2X  MZ . 30E " I Oqdp  . Aqdp 

2080  WAIT  2 

2085  GOTO  1970 

2090  ELSE 

2095  BEEP 

2100  PRINT  USING  "4X  ,HZ.3PE  ,2X ,MZ.30E":0qdp .Aqdp 

2105  Err=5P0 

2110  UA I T  2 

2  115  GOTO  1970 

2  120  ENO  IF 

2125  ENO  IF 

21301 

21351  LOOP  TO  SET  T,at 

2U0  IF    Ido-2    THEN 

2145  BEEP 

2150  INPUT  "ENTER  OESIREO  Tsal"  .Otld 

2155  PRINT  USING  "4X,""  PTsat   AT:st    Rate     Tv      Rat. 

2160  Oldl-0 

21G5  Old2=0 

2170  OUTPUT  709i"AR  AF33  AL35  UR5" 

2175  FOR  I-l  TO  3 

2180  Sun=0 

2185  OUTPUT  709i "AS  SA" 

2190  FOR  Jl-I  TO  2 

2195  ENTER  7P9.E1 iq 

2200  Sun»5um+El  in 

2205  NEXT  Ji 

2210  EUo-Sum/2 

2215  rid-FNTvsv(Elio) 

2220  IF  1  =  1  THEN  At  ld-Tld 

2225  JF  1-3  THEN  Tv-Tld 

2230  NEXT  I 

2235  IF  ABS(Atld-Otld >  -.2  THEN 

2240  IF  At ld'Otld  THEN 

2245  PfcFP  4000, .2 

2250  BEEP  4000,. 2 

2255  PEEP  4000,,  2 

2260  ELSE 

2265  BEEP  250.  .2 

2270  BEEP  250,. 2 

3275  PEEP  250  .2 

2280  ENO  IT 

2235  Errl-Al Id-Oldl 

2290  Oldl-Atld 

2295  Err2-Iv-01d2 

2300  01d2-Tv 

2305  PRINT  USING  "4X.5(MOO.DO  I  X  1  •  ;  Dt  1  d  A  t  I  d  Err  I  .Tv,  Err2 
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2310  WAIT  2 

2315  GOTO  21 70 

2320  ELSE 

2325  IF-  A8S<  At  Id  Ot  Id)  >.  I  THEN 

2330  IF  Atld'OHd  THEN 

2335  SEEP  3000,. 2 

2340  BEEP    3000,  .2 

2345  ELSE 

2350  BEEP   500,  .2 

2355  BrEP   800,  .2 

23G0  END    IF 

23G5  Errl-AUd-Oldl 

2370  Oldl-Alld 

2375  Err2=Tv-01d2 

2360  01d2-Tv 

2385  PPINT  USIN6  "4X  ,S< MOD . OD  ,  I  X >" i Dt 1 d  ,  At  1 d  ,Err  I  ,  Tv  ,Err2 

2390  UAIT  2 

2395  60T0  2170 

2400  ELSE 

2405  BEEP 

2410  Errl-At ld-Oldl 

2415  Oldl-Mld 

2420  Err2=Tv-01d2 

2425  ,0!d2=Tv 

2430  PRINT  USING  "4X  ,5(  MDO .  DO  ,  I  X  )"  iDt  Id  ,«Ud  ,Err  I  ,Tv  ,ErrZ 

2435  UAIT  2 

2440  GOTO  2170 

2445  ENO  IF 

2450  ENO  IF 

2455  END  IF 

2  450"  ERROR  TRAP  FOR  [do  OUT  OF  BOUNDS 

2465  IF  fdo>2  THEN 

2470  PEEP 

2475  GOTO  '890 

2490  ENO  IF 

24P5I 

24101  TAKE  DATA  IF  In-0  LOOP 

2495  PEEP 

2500  INPUT  "ENTER  BULK  OIL  V  Boo 

3505  OUTPUT  709i"AR  AF25  AL3E  VR5" 

2510  FOR  l^|  TO  12 

2SIS  OUTPUT  709i "AS  SA" 

2520  Sun-0 

252S  FOR  Ji-I  TO  2 

2530  ENTER  709i£ 

2535  Sun=Sum+E 

3540  NEXT  li 

2545  Enf (  I  >=Sum/2 

2550  NEXT  I 

2555  OUTPUT  709T'AR  AF62  At  63  "RS- 

2560  FOR  I=|  TO  2 

25G5  OUTPUT  709 i "AS  SA" 

2570  Sum  =  0 

2575  FOR  Il-I  10  2 

2580  ENTER  709 iE 

2585  Sum'Sum-tE 

2590  NEXT  J  i 

2595  IF    1=1     THEN   Vr-Sum/2 

2E00  IF    1=2    "HEN    tr-Sum/2 

2605  NEXT    I 

2610  ELSE 

2615  EflTEP    0Fi le2;Boo    Told*  ,Enf (  •  >  ,Ur  ,Ir 

2G20  END  IF 

2525' 

2E30I  COM'.'ERT  enf'S  TO  TEMP  "01.  T  .CURRENT 

2G35  rua»8 
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2640  FOR  1=1  TO  12 

2645  IF  Idtc.e  THEN 

:S50  IF  I-Ldtc 1  OR  I»Ldtc2  THEM 

2655  TU>— 99.99 

2SGO  GOTO  2710 

2665  ENO  IF 

2670  ENO  IF 

2E75  IF  1 1 1 • 4  THEN 

2680  IF  I>4  AND  I<9  THEN 

26P5  TUJ--99.99 

2E90  GOTO  2710 

2695  ENO  IF 

2700  ENO  IF 

2705  T( I  >-FNTvsv<Emf< I  I ) 

2710  NEXT  I 

2715  IF  1 1 1  <•  4  THEN 

2720  FOR  I»l  TO  4 

2725  IF  l=Ldtcl  OR  !=Ldtc2  THEN 

2730  Tuia-Twa 

2725  ELSE 

2740  Tuia  =  Tua+T(  I  ) 

2745  ENO    IF 

2750  NEXT    I 

2755  Tua-Twa/U-Idtc) 

2760  ELSE 

2765  FOR    I- I    TO    8 

2770  IF    I  =Ldtc I    OR    I-Ldt  c2    THEN 

2775  Twa=Tu,a 

2780  El  SE 

2785  Twa-Twa+TI  I  > 

2790  ENO    IF 

2795  NEXT    I 

2800  Twa-Tuia/<8-ldtc) 

2805  ENO    IF 

2810  Tld'f<9> 

2PI5  TldZ>T(l0) 

2920  TvT<  I  I  ) 

2925  IF    Itt    3    THEN 

2830  Tld2— 93.99 

2835  Tv-<  T(  I8HTI  I  I  )  1/2 

2940  ENO    [F 

2845  Tsnno'TI 12) 

2950  Anu-Ir 

2e55  Uolt»l'r«25 

2860  Q«Volt*Ano 

2eG5  IF  T 1 1  - 0  THEN 

2870  Kcu«FNKcu<Tua> 

2875  ELSE 

2e90  Kcu=Kcua(  Itn 

2985  ENO  IF 

28901 

2095'  FOI.irifh  CONDUCTION  E0UAT10N  UITH  CONTArT  RESISTANCE  NEGLECTED 

2900  Tu-Tua-0»L0e<D2/DI  >/(2«PI»Kcu«L  > 

2905  Thetab-Tw-I  Id 

2910  IT  rhetab  0  THEN 

2915  BEEP 

2920  INPUT  "TUALL  TSAT  [B-CONTINUE,  l=END>*,!e. 

2925  IF  Iev-0  THEN  GOTO  1875 

2930  IF  levl  IIIFN  3330 

29:5  END  IF 

29401 

29451  COMPUTE  "ARIOUS  PR0FFRTIE9 

2950  Tf  i  Jn-FMIf i |n( Tu.Tld) 

2955  Rho-FNPho(  If i In) 

2960  Mu-rtlMul  T  r  i  In  I 

2965  K»FNK( Tf i In ) 
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2370   Co=FNCo<  Tf 1  In) 

2375  Beta-FNBela(Tf ilw) 

2990  Ni-Mu/Rho 

2335   Aloha-K/(Rho«Co > 

2330   Pr-Ni/Alpha 

2335   Psat-FNPsat! Tld) 

3(il?0i 

30051  COMPUTE  NATURAL-CONUECTR'E  HEAT-TRANSFER  COEFFICIENT 

3010'  FOR  UNENHANCEO  ENO< S  ) 

3015   Hbar=l90 

30:0   Fe-(Hbar'P/(Kcu«A ) >" .S«Lu 

3025      Tanh-FNTanMFe  ) 

3030       Theta-Thetab'Tanh'Fe 

3035      X--<  9.81  'Ocla'Thetab-Oo    3«  Tanh/  (  Te«Ni  -A  Inha  )  )    .  I  SG5G7 

3040      Yy-<  l  +  <  .553/Pr  )"( 9M6 >)    (B/27) 

3045      Hbarc  =  K'/Do«(  .G  +  .387-X.-  'Yy  )~2 

3050   IF  ABS( (Hbar-Hbarc  >/Hbarc  P .001  THEN 

3055   Hbar-tHbar+Hbarc >».5 

30G0   GOTO  3020 

30S5   END  IF 

3070! 

30751  COMPUTE  HEAT  LOSS  RATE  THROUGH  UNENHANCED  END'S) 

3080   Q]=(Hbar«P»Kcu«A  >  . 5«Thetab« Tanh 

3085   Qc-O-2'Ql 

3030  As-PI»02*L 

3035' 

31O01  COMPUTE  ACTUAL  HEA1  FLUX  ANO  BOILING  COEFFICIENT 

3105   Qdn-Oc.'Ai 

3110  H<ube-OdD/ The tab 

31  IS' 

31301  RECORD  TIME  OF  OATA  TAKING 

3135   IF  lm*0  THEN 

3130   OUTPUT  709t"TD" 

3135   ENTER  709 1  To  Id* 

3140   ENO  IF 

31451 

31501  OUTPUT  DATA  TO  PRINTER 

3155   PRINTER  IS  701 

3IG0   IT  lov-a  IHFN 

3185   PRINT 

3170   PRINT  USING  ■  1 0X  ."-Data  Set  Number  -  " " ,000 ,2X , " "Bull  Oil  X    -  ""  ,DD. D  ,5X  ,  I 

4A' 1  J  Bod , TcldS 

3175   PRINT 

3180   PRINT  USING  "10X,""TC  No:     I       3       3       4       5       6       7 

8  - "  " 
3195  PRINT  USING  "  1 0X  , "" Tenp  :  ""  ,8<  IX  ,MOD.DD>"  )T<  I  >  ,T<  Z  >  ,T<  3  >  ,T<  4  >  ,T<5>  ,T(6  >  ,T< 
7  I  ,T(8  ) 

3190   PRINT  USING  "I0X,""  Tua     Iliad    Tllqd2   T.aor   P5al    Tjunp""" 
3135   PRINT  USING  " 1 PX  ,  2 < MOP . DO  .  I  X  )  ,  I  X  ,M0D . 00  .  I  X  ,2(1 X  ,M0D . 00 ) ,ZX ,MOD . D- ; Tua . T I d  , 
T  IdZ  ,Tv  ,Psat  Jsump 

3300   PRINT  USING  "I0X,""  The  tab   Htube      Qdo""" 

3305   PPINT  USING  "  I  0X  ,MPO.  30  ,  I  X  :MZ  .  30E  ,  I  X  ,MZ  .  3PE"  i  Thet  ah  Htube  ,Odo 
3210   END  IF 
3315   IF  Iov-I  THEN 
3220   IF  J-l  THEN 
3;35   PRINI 

3330   PPINT  USING  "10X  ""RUN  No  OilX   Tsat    Htube      Qdo        Thetab-" 
3235   ENO  IF 

3240   PRINT  USING  "I2X.3D  4X  DD  ,  3X  ,  MOD  .  DO  ,  3(  I  X  ,  MZ  .  3DE  ) "  1  J  .Boo  .  T  1  d  ,Ht  ube  ,QdD  ,  The  t 
ab 

3245   ENO  IF 
3250   IF  ln-0  THEN 
3255  BEEP 

32E0   IflFUl  "OK  TO  STORE  THIS  DATA  SET  (l=Y  0«N)7"  01 
32G5   ENO  IF 
3270   IF  01  -I  OR  In-  I  1HEN  J-.1H 
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3275   IF  Ok-I  AND  ln=0  THEN  OUTPUT  SFi le I 1B00 . TcldS  Enf<«):Ur,Ir 

3280   IF  In-I  OR  01=1  THEN  0U1FUT  9PlotiQdo  Thetab 

3285   IF  ln«0  THEN 

3230  BEEP 

3295   INFUT  "WILL  THERE  BE  ANOTHER  RUN  ( l-Y ,0-N  )1'  ,So_on 

3300  Nrun=J 

3305  IF  Bo.onOI  THEN  3330 

3310  IF  Go_on-I  THEN  Reoeat 

3315  ELSE 

3320  IF  J-.Nrun+l  THEN  Reoeat 

3325  END  IF 

3330  !F  In-0  THEN 

3335  BEEP 

3340  PRINT  USING  " 1 BX ," "NOTE :  "",ZZ.""  data  runs  uere  stored  in  file  " '  , I 0A ' i J - 

I  ,D2_file* 

3345  ASSIGN  ©r i 1 e I  TO  • 

3350  OUTPUT  9Fi le2iNrun-l 

3355  ASSIGN  SFilel  TO  D!_file* 

3360  ENTER  SFi  lei  tOate*  ,Ldtc  I  .LdtcC  ,1 1 1 

3365  OUirUT  9Fi IpriOatel  Ldlcl  ,Ldtc2  .1 1 1 

3370  FOR  1=1  TO  Nrun-I 

3375  ENTER  SFileliBop  Told*  ,Emf  (  •  >  ,Vr  ,Ir 

3300  Ol'TPUT  0Fi  le2iBop  ,Told*  Emf < • ) ,Vr , Ir 

3385  NE'<T  I 

3390  ASSIGN  9F  t 1-1  TO  • 

3395  PURGE  "DUMMi" 

3400  END  IF 

3405  BEEP 

3410  PRINT  USING  -|0X,"-NOTE:  "".ZZ,""  X-Y  oairs  uere  stored  in  plot  data  file 

""  I0A"U-I  P_fl!p$ 

54 15  ASSIGN  9Fi \r?    TO  • 

3420  ASSIGN  PPlot  TO  • 

3425  CALL  Stats 

3430  BEEP 

3435  INPUT  "LIKE  10  PLOT  DATA  (  I  -Y  ,9«N )?"  ,CH 

3440  IF  01  =1  THEN 

3445  CALL  Plot 

3460  EMP  IF 

3455  SUBf-NO 

34G0I 

3465 1  CUB'-'E  FITS  OF  PRUI'FRTY  FUNCTIONS 

3470  OFF  FNKcul  I  I 

3475'  OFHC  COFPER  250  TO  300  K 

34e0  Tl =1*273. 15      IC  TO  K 

34e5  K-4.34-.  1  1 2  •  T* 

34  90  RETURN  K 

3495  FNENO 

3500  OFF  FNMu<  T I 

35051  170  TO  360  K  CURUE  FIT  OF  UISCOUSI1V  . 

3510  TJ.  « 1  *273.  15      IC  TO  K 

3515  Mu=EXP< -4.4636+<  101  I  .  47'Tl  )  )•!  ,0E-3 

3520  RETURN  Mu 

3525  FNFND 

3530  OFF  FNCot  T) 

35351  IP0  TO  4P0  K  CURVE  FIT  OF  Co 

3510  Tl =T*273. 15      IC  TO  K 

3545  Cp>. 40186*1  .65007E-3M1  H  .5I494E-6-T1  ~2-6  .57RS3E- I 0» Tl     3 

3550  Co=Co«IP00 

3565  RETURN    Co 

3PG0  FNEND 

35G5  OFF  FNRhot  T  ) 

3570  Tl  -T+273.  15     !C  TP  1 

3575  X=l-<  I  .8«T1  /753.9S)   IK  ro  R 

35C0  Re- =  36.32*6  I  .  1464I4X  (  I  3  l  H  6  .  4  I  801  5  •  <  *l  7.  4.76  838  »X  "  .  SH  .  I  I  982  B-K  '2 

7585  Ro=R0' .062431 

3590  RETURN  Ro 
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3595      FNEND 

3EP0      DEF    FNPH  T  ) 

3605      Pr'FNCpl T J.FNMuf  T  >/FNK< T ) 

3610      RETURN    Pr 

36  I  5  FNEND 

3620      DEF    FNKl T) 

36251  T-360   l<    WITH    T    IN   C 

3G30  K-.  07I-. 000:61 »T 

3635  RETURN   K 

3640  FNEND 

3645  DEF    FNTanhiX  ) 

3650  P-EXP(X) 

3655  O-l/P 

3660  Tanh=<P-Q>/<P  +  0> 

3665  RFTURN    Tanh 

3570  FNEND 

3675  OFF    FNT.-3-.  <V) 

3E-F0  COM    iCc/   C(7)  ,Ical 

3585  T=L'iO> 

3690  FOR    1=1     TO    7 

3635  T=T+C<  I  l«V"I 

3700  NEXT    I 

37(15  IF    Ical-I     THEN 

3710  T-r-6.74?2934E-2<-T"<9.0377043E-3-T«<     9  .  33539  I  7E  -5  )  ) 

3715  ELSE 

3730  T=T+8.G2GB97E-2+TM3.7EI99E-3-T«5.0G89::S9E-5  I 

3735  END    IF 

3730  RETURN    T 

3735  FNEND 

3740  DEF    FN8eta( T) 

3745  Roo-FNRho(T* .  I  ) 

3750  Rom-FMRhoCT-.  I  ) 

3755  Beta-- 2/<Reo*Rom  >•<  Rou-Rom  )/ . 3 

3760  RETURN   Beta 

37G5  FNEND 

3770  DEF    FNTfllmt  Tw,Ild) 

3775  ffilm«<Tw»TldV2 

3790  RETURN    T"f  I  In 

37P5  FNEND 

3790  OFF  FUPsatl Tc  ) 

37951  0  TO  80  d<rg  F  CURVE  F I T  OF  Psat 

3600  Tf-I .B«Tr+33 

3805  Pa»5.94552S*Tf •<  .  I53520P2 H f •<  I  .  4840963E-3+ T f »9 . 6  I  5067 1 E-6  )  ) 

3ei0  Pg-Pa-14.7 

3815  IF  Pg  0  THEN       !  f-PSIG ,  —  in  Hq 

3G20  Psat'Pg 

3935  ELCE 

3F30  P5at=Pg-39.93.'l4.7 

3835  END  IF 

3840  PETURN  Psat 

3845  FNEND 

3?r^0  SU(3  Plot 

3855  COM  /CdJ//  fl(l0  10)  CM0)  B(  4  )  .Nop  ,  Iprnt  ,0do  ,  I  log 

3P60  INTEGER  Ii 

3865  PRIMTFR  IS  1 

3870  BFTP 

3875  INPUT     "I  IKE    DEFAULT    VA|  IJF5    FUR    PI  OT    (  1»Y  ,0-N>?"  ,Idv 

3030  eEFP 

3895  PRINT  USING  "4X  , " "Sc I «-r- 1  Option:"*" 

3890  PRINT  USING  "GX,"^  a    versus  delta-T 

3895  PRINT  USING  "SX,-"1   h  versus  de I tn  -T "  "  " 

3900  FRIMT  USING  "GX,""2   h  .prsus  a*"" 

3905  INPUT  On" 

3310  BFEP 

"915  INPUT  "SELECT  UNITS  ( 0«SI  ,  t -ENGLISH >"  tun 

3920  PRINTFR  IS  70S 
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3935  IF  Idv.t  THEN 

3930  BEEP 

3935  INPUT  "ENTER  NUMBER  OF  CYCLES  FOR  X-AXIS", C.< 

3940  BEEP 

3945  INPUT  "ENTER  NUMBER  OF  CYCLE?  FOR  Y-AXIS", Cy 

3950  PEEP 

3955  INPUT  "ENTER  MIN  X-VftLUF  (MULTIPLE  OF  l0)',Xnln 

39E0  BEEP 

3965  INPUT  "ENTER  MIN  Y-UALUE  (MULTIPLE  OF  l0)",Ynin 

3970  ELSE 

39/5  IF    Ocio=0    THEN 

39P0  Cy=3 

3965  Cx-3 

3?q0  <Min«. I 

3995  fnin-100 

4000  ENO  IF 

4005  IF  0oo=l  THEN 

4010  Cy*3 

4015  Cx»3 

4030  Xf,,n=.l 

4035  Ypum=I00 

40^0  ENO    IF 

4035  IF  Opo*2  THEN 

(040  rv=3 

4045  C-=3 

4P50  Xmin-1000 

4055  Ynin-100 

40F0  ENO  IF 

10F5  ENO  IF 

40-70  BFEP 

4075  PRINT  "INiSPIiIP  ^300. 3300  ,9300  .S900I " 

40P0  PRINT  "SC  0,100  0,l00iU  2,0i" 

40H5  Sf-100'C 

40Q0  SfytW/Cy 

4095  PRINT  "PU  0  ,0  PD" 

4  100  Nn-9 

4105  FOR  1=1  TO  r.  n 

4  110  Xat  Mniin  10  '  (  I-  M 

4  115  IF  I»C«+I  THEN  Nn> I 

4  130  FOR  J-l  TO  Nn 

4ln5  IF  1  =  1  THEN  PRINT  "  Tl  3  0" 

4130  IF  ]~2    THEN  FRTNT  "Tl   I  0" 

4  135  Xa-Xat».l 

4  I  4P  X»|.(!T(Xa/Xmin>»Sf< 

4  I  45  PRINT  "PA" iX  ,"  01  XT i " 

4  1F0  NEXT  J 

4  155  NEXT  I 

4IF0  PRINT  "PA  100  0iPUi " 

4IE5  PRINT  "PU  PA  0  0  PO" 

4  170  Nn-9 

4  175  FOR  1=1  TO  Cy+I 

4IP0  Yat-Ynin"  I0'(  T-l  ) 

4185  IF  [-Cy+I  THEN  Nn-I 

4IQ0  FOR  ).|  TO  Nn 

4195  IF  1  =  1  THEN  PRINT  "Tl  3  0" 

4300  IF  .1  =  7  IHEN  PRINT  "Tl   I  0" 

4305  Y*-Yat».l 

4310  Y»lRT(Ya/Ymln>-Sfy 

4315  FRINT  "PA  0  "|Y  "YT" 

4770  NFXT  J 

4325  NFXT  I 

4330  PRINT  "PA  0  ,  100  Tl  0  3" 

4335  Nn=9 

4  3  40  FOR  1  =  1  TO  f-  i  I 

4345  XafXmin-  10    t  t-l  > 

43r-0  IF     1=0.  4  1     THFN    Nn>  I 
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4255  FOR  J-l  TO  Nn 

42G0   IF  J-1  TMFN  PRINT  TL  0  2" 

42S5   IF  JNI  THEN  PRINT  "TL  0  I" 

4270   Xa-Xat»J 

4275  X«l  GT(  Xa/Xmln  )»Sf» 

4280  PRINT  "PA" iX ," , 1001  XT" 

42B5  NEXT  J 

4290  NEXT  I 

4295  FRINT  "PA  1 00  1 00  PU  PA  1 00  0  PD" 

4300  Nn-9 

4305  FOR  I-l  TO  Cv+I 

4310  Yat-Ynin«lf>"< I-l  ) 

4315   IF  I-Cy+I  THEN  Nn=l 

4320  FOR  J-l  TO  Nn 

4325  IF  J»l  THFN  PRINT  *TL  0  2 " 

4330   IF  J>l  THFN  FRINT  "TL  0  I" 

4335  Ya-Yat*J 

4340  Y-LGT(Ya/Ymin)»Sfy 

4345  PRINT  "PD  FA  I  00  , "  1  ,  " 1 T " 

4350  NEXT  J 

4355  NEXT  I 

4360  PRINT  "PA  I  00 .  I  00  PU" 

4.3(55  TFINI  "PA  0  -2  <:R  I  .5  2" 

4370  Ii-lGT(Xnin) 

4  375  FOR  1- 1  TO  C**  I 

4  380  Xa«Xmln« 10  <  1-1  ) 

4385  X=LGT(  Xa/Xnin  )«5f'. 

4310  PRINT  "PA"  iX  ,"  .01 " 

4395  IF  1 1 >=0  THEN  PRINT  "CP  -  2  ,-2 ;LB I 0iPR  -2  2  iLB"  i  I  i  :  " " 

44P0  IF  Ii  0  THFN  PRINT  TP  •  2  . -2  1 1.  B I  0 1  PR  B  ,2iLB"  1 1  i  I  *' 

4  405  I  1  =  1  i  H 

4410  NEXT  I 

4415  PRINT  "PU  PA  0  0" 

4420  I  i =L  PT( Ywm ) 

4425  Yie»ia 

4430  FOR  1=1  TO  Ty+I 

4435  Y*=Ymin»  10  (  I  -I  ) 

4440  Y«LGT(  Ya/Ynln)»Sfy 

4445  PRINT     "PA    0  ,  " | Y , " " 

44^0  PRINT    "CP    -4  ,-.2Sil.BI0;FR    -2 ,2 iLB" 1 1 1 1 " " 

4455  Ii-ti+1 

44G0  NEXT    I 

44G5  BFEP 

4470  INPUT    "WANT    USE    OFFAHIT    I  ABEt  9    (l«Y    0-NI?",IdI 

4475  IF    IrJI  (  .1     THFN 

4480  RrF.P 

4495  INPUT  "E.1IER  X -I. ABEL"  xiahelS 

4490  BFFP 

4495  INPUT  "ENTER  Y-LABEL"  ,Ylabel» 

4500  END  IF 

4505  IF  Ooo  2  THFN 

4510  PRINT  "5R  I,2iPUPA  40  -14;" 

4515  PRINT  "1R(T;PR  -1.9,3  fO  PR  1.2,0  PHiPP  .  5  ,  -4  ;  I.Buo  s  PR  .5,1  t " 

4530  PRINT  "LB-TiPR  .5  -  I  i  L  Rm  t  iFR  .S.li" 

4525  IF  Iun-0  THEN 

4530  PRINT  "LB  I  <K  >" 

4535  ELSE 

4540  FPINT  "IB)  (F)" 

4545  FUD  IF 

4550  FNP  IF 

4555  IF  Opo=2  THEM 

4';60  IE  [un-0  THFN 

45r.5  PRINT  "CR  I  .  5  ,  2  ;  PU  PA  40  .  l4ilBn  (U/ntSR  l,l.5iPR  0  .  5  .  I  il  R3  i  3R  l.5,2lPR 

0.5,-1  it  B)" 

4570  El SF 

4r,75  pnitu    "SR    I  .  5  , 2 1  PU   PA    34,-MilBn    (Btu/hnPR    .5  ,  .  S ;  I  B  .  :PR    .5, -.Si" 
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4?pfl  PRINT  "LPftiPR  .5,liSR  I  I.SiLBZiSR  I.S.ZiPR  .5-1 « L  B  >  i  ■ 

4595  END  IF 

4590  END  IF 

4595  IF  Ono  =  0  THEN 

4600  IF  lun  =  0  THEN 

4605  PRINT  "SR  I  .5  ,ZiPW  PA  -12  ,40iOI  0,  I  iLBci  (U/niPR  -I  ,0.5|5R  I ,! .5lLB2tSR  I 

,5,2iPR  I  ,.5iLB>" 

4GI0  ELSE 

4615  PRINT  "SR  I.S.ZlPU  PA  -IZ,3Z|DI  0,lil.8q  (Btu/hnPR  -  .  5  ,  .  5 :  LB  .  i  PR  .5, .Si" 

4620  PRINT  "LBftiSR  I,I.S|PR  -I,.5iLBZiPR  I  . . 5 1 SR  I.S,2iLfl>" 

4625  END  IF 

4630  END  IF 

4635  IF  Coo  0  THEM 

4640  IF     Iun-0    THFN 

4645  PRINT     "SR    I.S,2lPU    PA    -I2,38|DI    O.lilBh    (U/miPR    -l,.SlSR     I,1.5iLB2iSR     I. 

5,ZiPR    .5, .Si" 

4650  PRINT  "IB.iPR  .5  0|LBK>" 

4655  ELSE 

4650  PRINT  "SR  l.5,2lPU  FA  -IZ.ZBiDI  fl.lilBh  (Btu/hnPR  -.5  .  .  5iLB.  iPR  .5,.5i" 

4665  PRINT  "LPftiPR  -l,.5iSR  I.I.SiLBZlSR  I.S.ZlPR  .S,.5.LB.iPR  .5,.5;LBF>- 

4670  FNO  IF 

4675  END  IF 

4680  IF  Irll-'O  THFN 

4685  PRINT  "<">R  I.S.ZlPU  PA  50,-IG  CP~  i  -I  EN(  XI  abe  1  S  )  12  i  "0 1  LB "  i  <  1  abr  1  J  i  " " 

4690  PRINI  "PA  -14.50  CP  0  . " !-l EN< YlabelS >/Z«5/6i "01  0 , I  iLB"  1 1  label  $ i " " 

4695  PRINT  "CP  0.0  01" 

4700  FNO  IF 

4705  lnn«0 

4710  ReDeat ! ' 

4715  Kl  1-1 .E+G 

4720  Xul  — 1.E*S 

4725  len«0 

47-0  eEEP 

4735  INPUT  "UANT  TO  PLOT  DATA  TRON  A  FILE  ( I «Y ,0-N )'" ,01 

4740  IF  Ok-I  THEN 

4745  BEEP 

4750  IIIPIJT  "ENTEP  I  HE  NAME  OF  IMF.  DATA  FILE"  0_file* 

4755  ASSIGN  9F i 1 e  TO  D_f i 1 e* 

4760  BEEP 

47E5  PEEP 

4770  INPUT  "ENTER  THE  BEGINNING  RUN  NUMBER". Md 

4775  BEEP 

4780  INPUT  "ENTER  THE  NUMBFR  OF  X-Y  PAIRS  S ICRED" ,Npa irs 

4785  PEEP 

4790  INPUT    "CONNECT    DAIrt    WITH   LINE    <  I  »Y  ,0-N  >'"  ,  1c  1 

4735  BEEP 

4800  PRINTER  IS  I 

4805  FRINT  USING  "4X, ""Select  a  symbol:"*" 

4810  PRINT  USING  "6X,""I   Star     2   Plus  sign""" 

4815  FRINT  USING  "6X  "  "3   Circle    4   Square""" 

4820  PRINT  USING  "6X,""5   Rombus""" 

4825  PRINT    USING    "6X,""6      Rirjht-Slde-U0    triangle 

4830  PRINT    USING    "GX.""7      Uo-s  i  rle-doun    triangle""" 

4835  [MPUI  Svm 

4840  PRINTER  IS  705 

4845  FRINT  "PU  DP 

4OS0  IF  S/m-l  THEN  PRINT  "SM«" 

4855  IT  Sym-2  THEN  PRINT  "SMt" 

4660  IF  Svm  =  3  THCN  PRINT  "<":Mo" 

4865  IF  fid- 1  THEN 

4070  FOR  I-l  TO  (Md-I I 

4875  ENTER  £F i leiYa  Xa 

4  880  NT  XT  I 

4PP5  END  IF 

4890  FOR  t-1  Tn  Npalrs 

4895  ENTER  0F  i le i  *i  a  .Xa 
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4900  IF    Opo-1    THEN    Va=Ya/Xa 

4905  IF    Ono-2    THEN 

4910  Q=Va 

♦915  Ya-Ya/Xa 

4=520  Xa=Q 

4925  END    IF 

4930  IF    Xa<Xll    THEN    Xll»Xa 

4935  IF    Xa.>Xul    THEN    XuI>Xa 

4940  'IF    lun-l    THEN 

4945  IF    Ooo<2    THEN    Xa-Xa»I.B 

4950  IF    Ooo:0    THEN    Ya-Ya«.l7GI 

4955  IF    Odo-0    THEN   Ya=Ya».3l7 

49G0  IF    Odo-2    THEN    Xa-X*«.3l7 

4  9G5  END    IF 

4970  X-LGT(Xa/Xnin >«Sf • 

4975  Y=LGT( Ya/Ynin)«Sfy 

4980  IF  Sy«>3  THEN  PRINT  "SM" 

4985  IF  Syw  4  THEN  PRINT  "SR  I.  4, 2.4" 

4990  IF  Ic 1-0  THEN 

4995  PRINT  -PA"  X  .Y,"" 

5000  ELSE 

5005  PRINT  "PA"  .X,Y  ,"PD" 

5010  END  IF 

5015  IF  Syn  3  THEN  PRINT  "SR  1.2,1.6" 

5020  IF  Svn't    THEN  PRINT  "UC2  .  4  . 99  .0  ,  -8  , -4  ,0  ,0  , 8  4  ,0 . " 

5025  IF  Svm  =  5  THEN  FRINT  "UC3  0  .  99  ,  -  3  ,  -6  ,  -3  ,6  ,3  .6  .  3  ,  -5  i " 

5030  IF  Syn=G  THEN  PRINT  "UC0  .5 . 3  ,  99  .3  , -8  , -6  , 0  .  3  ,8 i " 

5035  IF  Svm-7  THEN  FRINT  "UCO  . -5. 3  ,99  ,-3  .8  ,fi  ,0  , -3  ,-8 i" 

5040  NEXT  I 

5045  PRINI  "PIT 

5050  BEEP 

5055  ASSIGN  SFile  TO  ■ 

50G0  XU-Xll/l  .2 

5055  Xul-Xul'1.2 

50701  GOTO  8040 

5075  END  IF 

5080  FRINT  "PU  5M" 

5085  ECET 

5090  INPUT  -WANT  TO  PLOT  A  POLYNOMIAL  < 1-Y  ,0«N >?" ,Go_on 

5095  IF  Po^on'l  THEN 

5100  8  FTP 

5105  PRINTER  IS  I 

5110  FRIMI  USING  "4X  . "  "Se  I  ec t  line  type:""" 

5115  PRIMf  USING  "6X,""0     Solid  line"'" 

5120  PRINT  USING  "SX,""I      Dasheii""" 

5125  PRINI  USING  "GX,""2,,  5  Longer  line  -  dash""" 

5130  INPUT  Ion 

5135  PRINTER  IS  705 

5140  RFEP 

5145  INFUT    "SELECT    (0-L.IN  ,  l-LOGI  X  ,Y  > )' ,Ileg 

5150  lDrnt  =  l 

5155  CALL  Pely 

5160  FOR  Xx-0  TO  Cx    STEP  C w 200 

5IR5  Xa=Xmin»10"X« 

5170  IF  Xa  XII  OR  Xa>Xul  THEN  5300 

5175  I.n=I-nt| 

5IP0  Pu-0 

5185  IT    Inn- I     THEN    tdf-Icn    MOO    2 

5190  IF    Inn=2    THEN    Irif-Ien   MOP-    4 

5195  IF     ton=3    THEN    [df-Icn    MOD    8 

5200  IF    Ipn-4    THEN    [df»Icn    MHO    IG 

5205  IF    tpn-5    THEN    [df-Icn    MUD    28 

52  10  IF    Irif  =  1     MIEN    Pu=l 

5215  IF    Ooe  =  0    THEN    Ya=f NPolyi Xa  ) 

5220  IF    Onc  =  2    AND    IIog'0    THEN    Y'a-Xa/FNPolv(  Xa ) 

5225  IF    Opo  =  2    AND    [  log-  1     fHFII    i  a-FNPol  v<  Xa  > 
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THEN 

Ya=Ya« 

.  1761 

THEN 

Ya' 

>Ya» 

.317 

THEN 

Xa- 

■  Xa« 

.317 

5230  IF    Ooo-I     THEN    Ya«FNPeJ v < Xa J'Xa 

5235  IF    Ya' Ymin    THEN    5300 

5240  IT    Iun-I     THEN 

5245  IF    Odo    2    THEN    xa=<a«l  .0 

5250  IF    Ooo   0 

5255  IF   Odo-0 

5260  IF    0po=2 

5265  END    IF 

5270  Y-LGTI Ya/Ymin)«Sfy 

5275  K-L6TCXa/Xi>Hn>»Sf« 

5290  IF    Y    0    THEN    r»0 

5295  IF  YM00  THEN  GOTO  5300 

5290  IF  Pu-0  THEN  PRINT  "FA" ,X , Y , "PO" 

5295  IF  Pu-I  THEN  PRINT  "PA"  ,X  .  Y  ,  "PU" 

5300  NEXT  X- 

5305  PRINT  "PU" 

5310  END  IF 

5315  eEEP 

5320  INPUT  "UANI  TO  QUI  1  <  I =Y ,0-N >"  ,  Iq t 

5325  IF  tqt»l  THEN  5335 

5330  GOTO  4715 

5335  PRINT  "PU  PA  0  0  SP0" 

5340  ciJBEND 

5345  OEF  FNHsmootMX  ,8oo  lsat  ) 

5350  DIM'  A<5>  B<5>  ,C(5)  .0(5  ) 

5355  DATA  . 20526  ,. 2S332  .. 3  I  9048  ,  .55322  .. 79909  1.00253 

5360  OATA  .74515  ,.72992 , .73199 , .71225 , .69472 . .64197 

5365  OATA  .  4  I  092  ,  . I  7735  . 25  I  42  ,  . 54806  .  . 8 1 91 6  ,  I . 0845 

5370  OATA  .71403  .73913  ,.73565,  .696691  . .665867, .61889 

5375  READ  A(  •  )  ,f)(  •  )  ,C(  •  )  ,D(  •  ) 

5380  IF    Boo    6    THEN    I=6op 

5385  IF    Boo-6    THEN    1-4 

5390  IF    Bop- 10    rilEN    1-5 

5395  IF    Isat-I     THEN 

5400  Hs»EXP(rtU  >+B<  I)«L06(X)) 

5405  ELSE 

5410  Hs-EXP(C(  I  >*DU  >»LOG<X>) 

54  15  END    IF 

5430  RETURN    Hs 

5425  FNENO 

54  !0  DEF  FNPolvfX  ) 

5435  COM  /Cply/  A(  I  0  ,  I  0  )  ,C(  I  0  >  ,B(  4  )  Noo  ,  Iornt  ,0po  ,  Hog 

5440  Xl*X 

5445  Pely-6<0> 

5450  FOR    I -I     TO    Noo 

5455  IF  llog-l  THEN  X1-L06IX) 

5460  Polv=Polv+B<  I  i'«l'I 

5465  NEXT  I 

5470  IF  [lcg-1  THEN  Pc  1  y =EXP< Po 1 y  ) 

5475  RETURN  Polv 

5480  FNENO 

5485  SUB  Poly 

5490  DIM  R<  10)  ,S<  10)  ,Sy<  12  >  ,S- i  12  >  ,x>.<  I00>  ,Yv<  100) 

5495  COM  'Colv/  A(  I  0  ,  I  0  )  CM  0  >  ,  B(  4  )  N  .  Iornt  ,0oo  .  1 1  eg 

5500  COM  /X-vy/  Xp(5),Yp<5) 

5505  FOR  1-0  TO  4 

5510  B(  I  )>0 

5515  NEXT  I 

5520  BEEP 

5525  INPUT  'SELECT  ( 0-FILE  I ■ KEYBOARD ,2=FR0GRAM  >", In 

5530  Im-rMI 

55  35  BEEP 

Sr-40  INPUT  "FN1ER  NUMBER  Of  x-r  PAIFS".No 

5545  IF  In-I  THEN 

5550  BEEP 

5555  INPUT  "ENTER  OATA  FILE  NAME"  0  file* 
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5SG0  BEEP 

5565  INPUT  "LIKE  TO  EXCLUDE  DA1A  PAIRS  U»Y  ,0»N  >?"  ,Ied 

5570  IF  led- I  THEN 

5575  BEEP 

5580  INPUT  "ENTER  NUMGER  OF  PAIRS  TO  BE  EXCLUDED"  .  Ipe» 

5585  END  IF 

5590  ASSIGN  PFile  TO  D_flle» 

5595  ELSE 

5S00  BEEP 

5605  INPUT  "WANT  TO  CREATE  A  DATA  FILE  ( I »V  ,0-N  )?'   ,  Yes 

5610  IF  Yes-I  THEN 

5615  BEEP 

SE20  INPUT  "GIUE  A  NAME  FOR  DATA  FILE"  ,D_file$ 

56^5  CREATE  BOAT  D_file*  5 

5630  ASSIGN  @File  TO  D_fileS 

5635  END  IF 

5640  END  IF 

5645  BEEP 

5650  INPUT  "ENTFR  THE  OROER  OF  POLYNOMIAL "  ,N 

5E55  FOR  1-0  TO  N»Z 

5660  5y(I)=a 

55C5  Sx<  I  1-0 

5670  NEXT  I 

5675  IF  led-'  AND  In-!  THEN 

5680  FOR  1  =  1  TO  Iue.< 

5685  ENTER  erUeiX  ,Y 

5690  NEXT  I 

5695  END  IF 

57PO  FOR  I- I  TO  No 

5705  IF  In- I  I  HEN 

5710  IF  Qpo-3  THEN  ENTER  8FlleiX,Y 

5715  IF  Opo'2  1IIFN  ENTER  PFileiY.X 

5730  IF  I  log-!  THEN 

5735  Xt-X'Y 

5730  X-LUGIXi 

5735  Y-LOG(Xt) 

5710  END  IF 

5745  END  IF 

5750  IF  tn-2  THEN 

5755  SEEP 

5760  INPUT  'ENTER  NEXT  X-Y  PAIR",X,Y 

5765  IF  Yes»l  IHEN  OUTPUT  9FiletX,Y 

5770  END  IF 

5775  IF  In  3  THEN 

5790  X-(  I  )-X 

5785  Yv(I)-Y 

5790  ELSE 

5795  X  =  Xp<  I- I  ) 

5800  Y  =  Yp<  I-l  ) 

58P5  END  IF 

58I0  RI0J-Y 

58I5  Sy(0  )  =  Sv(0  )H 

5830  S(  I  )-X 

5835  S-  i  I  )«Sx(  I  )+X 

5830  FOR  1=1  TO  N 

5835  R<  J  )-R( J-l  )»X 

5840  S»( J  )  =  Sy(  J  >  +  Ri  J  I 

5845  NEXT  J 

5850  FOR  J-3  TO  N»2 

5855  S(  J  >  =  S<  J-  I  )»X 

5860  5.-  i  J  »-Sx(  J  )iS<  .!  ) 

5065  NEXT  J 

5870  NEXT  I 

5875  IF  Yes-I  AND  In-3  THEN 

5880  BEEP 

5035  PRINT  USING  "  I2X  ,00 


X-Y  oaln  were  stored  in  file 


I0A" iNp  ,D  file* 


103 


5890 
5895 
5900 
5905 
5910 
5915 
5970 
5915 
5930 
5935 
5940 
5345 
5950 
5955 
59S0 
53G5 
5970 
5975 
5980 
5995 
5990 
5995 
6000 
6005 
SO  I  0 
6015 

gc:0 

G025 

6030 

G335 

EPANC 

6040 

S045 

6050 

SP55 

6060 

G065 

6370 

6075 

5090 

G095 

G090 

6095 

6100 

6  I  PS 

Rl  10 

61  15 

6  i:0 

6  I  35 

6  130 

6125 

G!  40 

GI45 

CIS0 

6iS5 

6160 

6155 

GI70 

EI75 

6190 

GI8S 

6190 

GI9S 

6300 

G305 

6210 


END  IF 

S-.  (0  >-No 

FOR  1=0  TO  N 

C(  1  >-Sy<  I  ) 

FOR  J=0  TO  N 

A<  I  ,J  >-S*<  I  Hj  ) 

NEXT  J 

NEXT  I 

FOR  I -8  TO  N-l 

CALL  Dividei  I  ) 

CALL  Subtract! I  +  l  ) 

NFXT  I 

B(M)-C(N  >/A<  N,N  ) 

FOR  l. a    TO  N-l 

eiN- 1  -1  >-c<n-i-i  ) 

FOR  J-0  TO  I 

B<N-I-I>»8<N-I-I  )-A(N-l-I,N  J  )>B(N-J  ) 

NEXT  J 

0<  N-  I  -I  )-B(N-  I  -I  )/A(N-l-l  ,M-  I  -I  ) 

NEXT  r 
'PRINTER  IS  701 
IFF  I  NT  B(  •  ) 
'PRINTER  IS  705 

IF  lornt*0  THEN 
I  2X 


EXPONENT     COEFFICIENT"" 
I5X  ,DD  ,5X  .ND.70E"  I  I  ,B(  I  ) 


PRINT  USING 
FOR  1-0  TO  N 
PRINT  USING 
NEXT  I 
PRINT  "  " 

PRINT  USING  "I2X.""DA1A  POINT       X 
'"  "  • 

FOR  r«l  TO  No 
Yc-B(8> 
FOR  J -I  TO  N 
>c-Yc»B<  J  >-X-.<  I  >  'J 
NEXT  J 
D"Yvl  I  )  Yc 

PRINI  USING  "  ISX  ,3D.4X  . 4  (  MD . 5DE  .  IX  ) "  ;  I  ,Xx<  I  )  ,Y> 
NET  I 
END  IF 

ASSIGN  8F1 le  TO  • 
SUBfTNO 

SUB  Divide! M  1 
COM  'Colv'  A<  10  . 
FOR  I-M  TO  N 
Ao-A(  I  ,M  ) 
FOR  J=M  TO  N 
A(  I  J  )--A(  I  .J  >/Aa 
NEXT  J 

C<n-C(  I  )-'Ao 
NEXT  I 
5U8EN0 

Sl'B  Subtract  It  > 
COM  /Coly/  A!  10, 
FOR  I-K  TO  N 

J-K-l  TO  N 

J  >«A(K-I  J  )-A(  I 
J 

I  )  Cl  I  > 


YICALCULATED  )  DISCR 


(  I  J.Yc.O 


0  )  ,C(  10  >  ,0(4  >  ,N  ,Iprnt  ,0po  ,  I  log 


10).C(I0)    0< 4 >  ,N . Iprnt  ,0do  ,  I  leg 


FOR 

A!  I    J 

NEXT 

C(  I  l-C<K 

NEXT     I 

SU6FN0 

SUB    Flln 

CUM    'Col//    A!  10,  10  )  ,C<  10 )    B( 4  > 

COM    /X.  ,v/     <n(5>    Y»<5) 

PRINTER  IS  705 

BEEP 


N  ,  Iornt  ,0po  -  I  log 
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621  5   INPUT  "SELECT  <0»h/h07.  same  tube  ,  I -h<  HF  >/h<  sm  ) "  ,  Irt 

6220  beep 

6225  input  "uhich  T?at  ( 1-6.7,0— 2. 2>"  ,isai 

6230   Xnun  =  0 

G23S  Xma.-.-IO 

6:i0  Xsteo-2 

G245   IF  Irt-0  THEN 

6250  Ymln-0 

6255  Yna*-I .4 

62G0  Ystep-,2 

r.265  ELSE 

6270  Ynm  =  0 

G275  Ynax-15 

62e0  Ystep-S 

6295  END  IF 

6290  BEFP 

6295  PRINT  "  IN i SPl ; IP  2300  ,2200  ,8300  .6900 i" 

6300  PRINT  "SC  0.100,0  I00;TL  2  0 i" 

6305  Sfx"l00/(Xnax-Xnin) 

G3I0  Sfy-l00/< Ymax-Ynln) 

6315  PRINT  "PU  0  0  P0' 

6320  FOR  Xa-Xnin  TO  Xnax  STEP  Xsteo 

6325  X»( Xa-Xnin )»Sfx 

6330  PRINT  "PA"  ;X  ,"  .01  XT i" 

G335  NEXT  Xa 

£340  PRINT  "FA  !00,0iPUi" 

6345  PRINT  "PU  PA  0  0  PO" 

6350  FOR  Ya'Ymin  10  rna-  STEP  Ysteo 

6355  Y  =  <  Va-Ynm  >»Sfy 

63G0  PRINT  "FA  0  ." ;Y  ,"YT" 

6365  NEXT  Yj 

G370  PRINT  "PA  0  ,  1 00  TL  0  2" 

6375  FOR  Xa=Xmin  10  Xroa-  STEP  Xstep 

S380  X-< Xa-Xnin )«Sf a 

6395  PRINT  "FA" iX,"  ,100!  XT" 

G390  NEXT  Xa 

6395  PRINT  "PA  100,100  PU  PA  190,0  PO" 

E40O  FOR  Ya-Ynin  TO  Ynax  STEP  Ysteo 

6405  Y=<  Ya-Ynln>»Sfy 

6410  FRINT  "PO  FA  I  00  ,  "  .  Y  , " Y T " 

6415  NEXT  (a 

6420  PRINT  "PA  I  00  ,  I  00  PU" 

6425  PRINT  'TA  0.-2  SR  1.5,2" 

6430  FOR  Xa-Xmin  10  Xna.,  STEP  Xsteo 

6435  X'<Xa-Xnin)»Sf- 

5  4  4  3  FRINT  "PA" ; X . " ,0 : * 

6445  PRINT  "CP  -2  ,-|  iLB" ;Xai"" 

6450  NEXT  Xa 

6455  PRINI  "PU  PA  0,0" 

6460  FOR  ra=Ymin  10  Ynax  SIEP  Ysteo 

6465  IF  ABS(Ya)  I.E-S  THEN  Ya-0 

5470  Y  =  ( Ya-tnln  >»Sfy 

6475  PRINT    "PA    0  ,  " ; Y , " " 

6480  PRINT  "CP  -4  ,-  .  25iLB"  I  Ya :  " " 

6455  NEXT  1a 

6490  Xlabell-"Oil  Percent" 

6495  IF  Irt-0  THEN 

E5O0  YlabelS-"b/hOV 

6505  ELSE 

6510  y  labelS»"h/h*nooth" 

6515  END  I^ 

5520  PRINT  "SR  l.5,2iPU  PA  S0  -10  CP"  ;  -L  £N<  X  I  abel  ?  >  '3  i  "9  :t  B"  :  X  label  J  i" 

6525  PRINT  "PA  -11,50  CP    0  ," i-LEN< Ylabel$  )/2«5/Si "Dl  0  , 1 : 1 B"  i  Y 1 abel 1 1  ' 

G530  PRINT  "CP  0  0" 

6535  Ipn-0 

R540  BFFP 
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6545  INPUT  -UANT  TO  PLOT  DATA  FROM  A  FILE  (  I -Y  .0=N  )T  ,Q\  o 

6550       lcn*0 

5555  IF    01  o-l    THEN 

6560  BEEP 

G565  INPUT  "ENTER  THE  NAME  OF  THE  DATA  F ILE '  ,D_f I  1 eS 

6570  BEEP 

6575  INPUT  "SELECT  (0-LINEAR.  1  =LOG< *  ,Y >"  ,  I  log 

6580  ASSIGN  SFile  TO  D_file* 

6585  BEEP 

6590  INPUT  "ENTER  THE  BEGINNING  RUN  NUMBER"  ,Md 

6535  BEEP 

6600  INPUT  "ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED-  ,N0alrs 

66P5  BEEP 

6610  INPUT  "ENTER  DESIRED  HEAT  FLUX-,0 

6615  BEEP 

6620  PRINTER  IS  I 

S625  PRINT  USING  -4X, "Select  a    synod:""" 

6630  PRINT  USING  "4X,""1   Star    2   Plus  sign"- 

6635  PRINT  USING  "4X.--3   Circle  4   Souare""" 

6640  PRINT  USING  "4X,"-5   Rombus"-" 

6645  PRINT  USING  "4X,--G   R  I  ght  -  s  i  de-uD  tnangle--- 

E650  PRINT  USING  "4X,"-7   Uo-slde-doun  triangle-"" 

6655  INFUT  Syn 

6660  PRINTER  IS  705 

6665  PRINT  'PU    OI- 

6670  IF  Svm=l  THEN  PRINT  "SM-" 

6675  IF  Svm=2  THEN  PRINT  "SN+" 

6GeO  IF  S>m=3  THEN  PRINT  "SMo- 

6685  Nn-4 

6690  IF  I  log-l  THEN  Nn^- 1 

6695  IF  Md.'  1  THEN 

6700  FOR  I-l  TO  (Md-I  > 

6705  ENTER  ?Fl ieiXa ,Ya 

5710  NEXT  I 

S715  ENO  IF 

6720  Q1»Q 

6725  IF  I log-l  T MEM  0-LOG(Q  ) 

6720  TOR  I-l  TO  Noairs 

6735  ENTER  3File;Xa  B(  • ) 

6740  Ya  =  R(0) 

6745  FOR  K-l  TO  Nn 

6750  Ya-Ya+B(K  1-0'K 

6  755  NEXT  K 

6~G0  IF  I  log- 1  THEN  Ya-EXP(Ya) 

6765  IF  Ilog-0  THEN  Ya=01/Yrt 

6770  IF  Irt-0  THEN 

6775  IF  Xa-0  THEN 

6780  Yo=Ya 

6795  Ya-I 

6790  ELSE 

6795  Ya»Ya/Yo 

6800  END  IF 

6805  ELSE 

6910  Hsn-FNHsmootM  Q  ,Xa  ,  Isat  > 

6815  Ya-Ya/Hsm 

6B20  ENO  IF 

6825  X.<<  I-l  >-Xa 

6830  Yy< I-l  )-Ya 

6835  X"<Xa-Xmtn  >«Sf* 

G840  V<  Ya-Yminl-Sfv 

6845  IF  Svrt  3  THEN  PRINT  "SH" 

G050  IF  Svn  4  TIICN  PRINT  "SB  1.4,2.4" 

6655  PRINT  "PA"  X  >   " 

6860  IF  Svn  3  THEN  PRINT  "5R  1.3,1.6" 

6665  IF  Sym  =  4  THEN  PRINT  "UC2.4.99  0,-8,-4,0,0,8.4,0,1" 

6870  IF  Svn=5  1  HEN  PRINT  "UC3  .0  99  . -3  . -6  . -3  ,6  ,  3  .6  .  3  . -6  i  " 
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6875  IF  Svn>8  THEN  PRINT  "UC0 ,S . 3 , 99  ,3  ,-8  ,-6  .0  ,3  8  i" 

G880  IF  Svn-7  THEN  PRINT  "UC0 , -5 . 3 . 99  ,-3  , 8  .6  ,0  , -3  . -8 < ' 

S885  NEXT  I 

6093  BEEP 

6895  ASSIGN  (File  TO  • 

6  900  END  IF 

6905  PRINT  "PU  SM" 

6910  BEEP 

6915  INPUT  "UANT  TO  PLOT  A  POLYNOMIAL  C 1 «Y  ,0-N  >7"  ,0V p 

6920  IF  01 o-l  THEN 

6925  BEEP 

6930  INPUT  "SELECT  ( 0-L INEAR  .  I =L06( X  ,  Y  )  )"  ,  I  log 

6935  Iprnt-I 

6940  CALL    Poly 

6945  FOR    Xa»Xmn    TO    Xna*    STEP    Xsteo/25 

6950  Icn-Ien+1 

6955  Ya»FNPoly(Xa) 

E950  Y-(  Ya-Ynin  >«Sfy 

6965  X  =  (  Xa-Xnin  i-Sfx 

6970  IF    Y    0    THEN    Y«0 

6975  IT  Y  100  THEN  GOTO  70:5 

6980  Pu=0 

E99G  IF  Ion»l  THEN  [df»Icn  NOD  2 

6990  IF  Ipn-2  THEN  Idf=Icn  MOD  4 

6995  IF  Ion  =  3  THEN  Idf-[cn  MOD  8 

7000  IF  Ion-4  THEN  Idf-Icn  MOO  16 

7005  IF  Ion=-5  THEN  Idf-Icn  MOO  32 

7010  IF  Idf»l  THEN  Pu  =  l 

7015  IF  Pu=0  THEN  FR1NT  "PA"  ,X  1    "PD" 

7020  IF  Pu«l  THEN  PRINT  "PA"  X,Y,"PU" 

7025  NEXT  Xa 

7030  PRINT  "PU" 

7035  Ipn«Ipn+1 

7040  GOTO  6540 

7045  END  IF 

7C50  BEEP 

7055  INPUT  "U'ANT  TO  QUIT  (  I «Y  B-N  }?"  ,  Iqui t 

7060  IF  tquit-l  THEN  7070 

7065  GOTO  6540 

7070  PRINT  'PU  SP0" 

7075  SUGPND 

7080  SUB  Stats 

7085  PRINTER  IS  701 

7090  J-0 

7095  K=0 

7100  BEEP 

7105  INPUT  -PLOT  FILE  TO  ANALYZE?"  , Ft le* 

7!  10  ASSIGN  @Fi  le  TO  I-!  leS 

7115  BEEP 

7120  INFUT  "LAST  RUN  No M 0-OUI T  ) "  .Nn 

7125  IF  Nn=0  THEN  7305 

7130  Nn-Mn-J 

7135  Sv-0 

7140  Sv=0 

7145  Si'0 

7150  Sxs-C 

7155  Sys-0 

7IR0  Srs=0 

71S5  FOR  1=1  TO  Nn 

7170  J-JH 

7175  ENTER  VFileiO.T 

7180  H-O/T 

7185  S.--S-+0 

7190  Sxs«S<s+Q  2 

7195  Sy-Sy+T 

7200  Sv5-Svs<-T  2 
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7205  Sz-Sz+H 

7210  Szs-Sza»H-2 

7215  NEXT  I 

7220  Qave«S*/Nn 

7225  Tave-Sy/Nn 

7230  Have-Sz/Nn 

7235  Sdevq-SQR( ABS<  (  Nn-S-  S-Sx  2  »/  (  NnX  Nn-  I  )  )  )  ) 

72  40  3devt«5QR< A8S<  (Nn-Svs-Sy  2  >/(Nn«(Nn-l  )  )  )  ) 

7245  Sdevh-SQR< ABSf (Nn'Srs-Sz  2  ) ' < Nn«< Nn- I  ) ) ) ) 

7250  Sh=l00«9devh. Have 

7255  So-100'SdevQ/Qave 

7260  St*l0O«Sde*t/Tav.e 

72B5  IF  K-l  THEN  7295 

7270  PRINT 

7275  PRINT  USING  ■ 1  I  X  , " "OATrt  FILE: " "  ,  I4A" iFi le* 

72S0  PRINT 

7285  PRINT  USING  MIX, ""RUN  Htube     SdevH   OdQ        SdevO    Thetab  SdevT" 

7290  K-l 

7295  FRINT  USING  "  U  X  ,DO  ,  2<  2X  ,D  .  3DE  ,  I  X  ,  3D  .  2D  >  .  2X  .DP  .  3D  ,  I  X  ,3D  .  2D"  ;  J  .He /- ,5h  ,Oave 
,Sg ,Tave  ,  S  t 

7300  GOTO  7115 

73P5  ASSIGN  9Filel  TO  • 

7310  PRINTER  IS  I 

7315  SUBENO 

7320  SUB    Coef 

7325  COM  /Cply/  AU0.I0)  CM  0  »  ,B(  4  >  .N  .  Ipmt  ,0oo  .  Ilog 

7330  PCEP 

7335  INPUT  "GIVE  A  NAME  FOR  CROSS-PLOT  FILE-,Cpf* 

7340  CREATE  BOAT  Cof*  ,2 

7345  ASSIGN  9F 1 le  TO  Co  f5 

7350  BL'EP 

7355  INPUT  "SFLECT  IB-LINEAR  ,  I -I  0G<  X  ,Y  >>"  ,Uog 

73B0  BEEP 

73B5  INPUT  "ENIER  OIL  PERCENT  (  - 1 -STOP  I "  ,Boo 

7370  IF  Bop  0  THEN  7390 

7375  CALL  Poly 

73B0  OUTPUT  9Fi leiBop  ,B<  •  ) 

7395  GOTO  7360 

7390  ASSIGN  SFile  TO  • 

7395  SUBENO 
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APPEHDII    C 
COHPOTEB-CONTHOLIED    VALVE   PEOGEAM 

The  computer- control led  valve  was  controlled  by  an 
Octagon  Systems  SYS-2A  microcontroller.  The  SYS-2A  is  a 
complete  computer  system  on  a  single  card  and  requires  only 
a  5  volt  supply  for  operation-  The  SYS-2A  has  provisions  for 
4,  2-5  volt  A/D  data  inputs  and  8  high  current  digital 
outputs.    The  card  was  configured   as   shwon   in  Table    4    below. 

TABLE    4 
Configuration  of   SYS-2A   I/O  Channels 

Data   Input   Channel  Purpose 

3#AiO   =    8  Valve   position 

3#A10   =   9  Liquid   R-  1 14    Temperature 

Data   Output    Channel  Purpose 

3#A00  Shut   Valve  Relay    (0=stop, 1=shut) 

3#A01  Open   Valve   Relay    (0=stop,  1=open) 

The  proportional-integral-derivative  control  program 
below  was  written  in  NSC  "tiny  BASIC."  The  SYS-2A  microcont- 
roller operator's  manual  includes  an  appendix  of  tiny  EASIC 
commands-  Remarks  are  provided  below  that  were  not  in  the 
actual  program  (to  minimize  execution  time)  for  clarity.  The 
program  was  stored  in  EEPROM  added  to  the  card  after  program 
debugging   was   complete. 

10    Let   A=15:Let    B=5: Let  C=10:Let   F=0:Let   1=0      Initialize 

15    Let    5)#A00=0:Let    a)#A01  =  0  Stop   Valve 

20    Print  "Set   Constants  (1  =  Y,  0  =  N)  ?":Input    Q 

30    If  (Q=0)then    GOTO    60 

40    Print"A,B,C,G?":Input    H,J,K,L 

50   Let    A=H:Let    B=J:Let  C=K:Let   G=L 

60    Prinf'Input   Desired  Temp (-1=End) ":Input    R  Enter    Temp 
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70    Let    o)#A  10=9: Delay    2:Let    M=1  6*3#A  1 1+3# A1  2/16      Read    Channel    9 

80    Let    o)#A10=8:Delay    2:Let   P=16*S#A 1 1+3#A 12/16      Read   Channel    8 

90   Print"Temp="/M,"Pos=",P 

100    If  (R<0) then    GOTO    800  Exit   Program 

200    Let    o)#A10=8:Delay   2:Let    P=16*a)#A  1 1+o)#A  12/1  6    Read   Channel    8 

210    Let    S=0  Ave    10    Readings 

220    For    N=1    to    10    Step    1 

230    Let    a#A10=9:Delay  2:Let    H=16*a#A1U3#Al2/1 6   Read  Temp 

240   Let   E=M-R:Let   S=S+E  Error; average 

250    Next   N 

260    Let   E=S:Let    D=S-F:Let   F=E  Error; Derivative 

262    If  (E<0)then   1=1-1  Integral 

264    If(E>0)then   1=1+1  Integral 

266    If  (K-31000)  then    I=-31000  Binary   Limit 

268    If (I>31000)  then    1=31000  Binary    Limit 

270    Print"M=",M,  "E=",F,  "D='\  D,  »I=",I 

300    If  (E<-1000) then    GOTO    600  Control    Band 

310    If (E>2000)  then    GOTO   700  Control   Band 

400    Let    V=  (E/A)+  (B*D)  +  (I/C)  PID    Value 

410    Print"V=", V,"ET="/E/A,"DT=",B*D,"IT=",I/C         Debug    Tool 

420    If(V>0)then    GOTO    500  Above    Desired 

430    Let    V=V*(-1)*G  Shut    Faster 

440    If(V<1)then    V=  1  Check 

450    If  (V>1040)  then    7=0        Delay    Function    limited    1-1040    entry 

460    If  (P<150) then   V=15        Slow   Valve    near   end    of    travel 

465    If(P<50)then    V=10 

470    If{P<20)then    GOTO  200      End   of   Travel   Limit 

480    Let    a)#A00=1:Delay  (V)  :Let    S#A00=0  Shut   Variable   Time 

490    GOTO    200 

500    If  (P<2500) then   GOTO    200      End    of    Travel   Limit 

510    If(V<1)then    V=1  Check 

520    If  (V>1040)  then    V=0        Delay    Function    limited    1-1040    entry 

530    Let    a)#A01=1:  Delay  (V)  :  Let    3#A0  1  =  0  Open    Variable    Time 

540    GOTO    200 

600    If  (P<20)then    GOTO  200      End   of    Travel   Limit 
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610  Print"600  Loop"  Out    of   Bound    Warning 

620  Let    V=0  Max   valve   shutting  speed 

630  If  (P<250) then   V  =  50  Slow  near   end    of    travel   limit 

640  If  (P<150)  then    V=15 

645  If(P<50)then    V=10 

650  Let    5)#A00=1:  Delay  (V)  :  Let    o)#A00=0  Shut    Valve    Fully 

660  Let    5)#A10=8:  Delay  2:Let    P=  16*o)#A1 1 +3#A12/1  6   Read   Channel    8 

670  GOTO    600 

700  If (P>2500)  then    GOTO    200      End   of    Travel   Limit 

710  Print"700   Loop"  Out    of   Bound   Warning 

720  Let    5>#A01=1:  Delay   0:Let    o)#A01  =  1  Open   Valve   Fully 

730  Let    3#A10=8:Delay  2:Let    P= 16*3#A1  H-S#A12/1 6   Read  Channel    8 

740  GOTO    700 

800  END 
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APPENDIX  D 
EXAHPLES  OP  REPRESENTATIVE  DATA  EONS 

The   two     printouts  below      are  samples      of   the      output  of 

DRP2.         The   first    printout    (HH79)       is  for    the    smooth    tube   at 

high    heat      flux    with      0   percent     oil.        The      second   printout 

(HF101)    is   for   the   High   Flux   tube   under   similiar    conditions. 


Month,  date  and  tine  :03: 2S i I  0:02 : I  8 

NOTE:  Program  name  :  DRP2 
Disk  number   ■  02 
Old  file  name:  WH79 

This  data  set  taken  on  :  02 : I  0: I  2 :S0: I 5 
Tube  Tvoe:      4 


Data  Set  Number 


I   Bulk  Oil  X  -  0.0     02: 10: 13:16:18 


TC  No:     I       2       3       4       5       S       7       3 
Temo  :   IG.42   14.69   14.61   15.00   14.67   14.91   15.79   14.74 

Tua     Tliad    Tliad2   Tvaor   Psat    Tsumo 

IS. 23   -2. 25    -2.20    -.63   -6.31   -17.7 

Thetab   Htube       Qdo 

1G.935   S.704E+-03   9.660E+04 


Data  Set  Number 


2   Bulk  Oil  J  -   0.0      02: 10: 13: 19:35 


TC  No:     I       2      3      4      5      6      7      3 
Temo  :   16.49   14.60   14.54   14.93   14.52   14.74   16.69   14.57 

Tua     Tliad    Tliod2   Tvaor   Psat    Tsumo 

15.15   -2.27    -2.10    -1.44   -6.32   -17.1 

Thetab   Htube       Qdo 

16.866   S.74SE+03   9.689E+04 

Data  Set  Number  -    3   Bulk  Oil  !  •   0.0      02:10:13:19:56 

TC  No:     I       2       3       4       5       6       7       8 
Temo  :   16.48   14.56   14.52   14.95   14.50   14.57   16.62   14.67 

Tua     Tliad    Tliad2   Tvaor   Psat    Tsumo 

15.12   -2.27    -2.13    -1.51   -6.32   -17.0 

Thetab   Htuoe       Odo 

16.840   S.760E+03   9.701E+04 
NOTE:  03  X-Y  oairs  were  stored  in  plot  data  file  PUH79 

DATn  FILE:PUH79 

RUN  Htube      SdevH    Qdo         SdevQ    Thetab  SdevT 
3   5.736E+03     .51   9.683E+04     .22   16.980     .29 
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Month,  date  and  tine  : 03 : 25 : t 0 : 03 : 26 

NOTE:  Program  name  :  DRP2 
Ox  5k  number   ■  02 
Old  file  name:  HF 1 0  I 

This  data  set  taken  en  :  02:16:19:06:52 
Tube  Type:      5 

Data  Set  Number  -    I   Bull:  Oil  Z    >   0.0     02:16:19:31:02 

TC  No:  I  2       3  4      5      6      7      8 

Temo  :  5.52  2.63    3.76  4.57  3.93    4.52    7.31    4.17 

Tua  Tliqd  Tliqd2   Tvapr  Psat    Tsumo 

4.55  -2.14  -1.87      .37  -6.20   -15.7 
Thetab  Htube      OdD 

3.370   2.80IE+04   9.440E+04 

Oata  Set  Number  -    2   Bulk  Oil  X  -   0.0     02:16:19:31:16 

TC  No:     1  2      3  4.5      6      7      8 

Temo  :    5.54  2.65    3.76  4.57    3.93    4.53    7.31    4.18 

Twa     Tliqd  Tliqd2   Tvapr  Psat    Tsump 

4.56  -2.14  -1.90     .32  -6.20   -15.6 
Thetab   Htube      Odp 

3.379  2.792E+04   9.433E+04 

Data  Set  Number  -    3   Bulk  Oil  X    -   0.0      02:16:19:31:27 

TC  No:     I       2      3      4      5      6      7      8 

Temo  :  5.55    2.63    3.76    4.57    3.93    4.53    7.31    4.17 

Tua  Tliod    Tliod2   Tvaor   Psat    Tsumo 

4.56  -2.14    -1.91      .26   -6.20   -15.5 

Thetab  Htube       Odo 

3.369  2.807E+04   9.4S6E+04 
NOTE:  03  X-Y  pairs  were  stored  in  plot  data  file  PHF10I 

DATA  FILE:PHFI01 

RUN  Htube      SdevH    Odo        SdevO    Thetab  SdevT 
3   2.800E*04     .27   9.443E+04     .12    3.373     .16 
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APPENDIX    E 

ONCERTAINTY    ANALYSIS 

The  uncertainty  of  the  heat-transfer  coefficient  at  37 
kW/m*  and  5  kW/m2  of  runs  WH79  and  HF10  1  are  analyzed  below. 
The  analysis  is  based  on  the  Kline-McClintock  [Bef.  24] 
method   of    uncertainty  analysis. 

The   heat-transfer  coefficient   is: 


h  = 


(E.1) 


"wo        "^sat 


and 


T 


-T  =  Twj    -   Qc    in    (D2/Dl)        _T  (E.2) 


wo  sat  2   it  k   L  sat 

where 

h   =    heat-transfer    coefficient 
q      -    heat    flux    corrected   for    end    losses 
Two    =   average  outer    wall   temperature 
T  sa-t-    =    saturation   temperature 
^wi    =    avera9e  inner    wall   temperature 
Qc   =   heat    input   corrected   for    end   losses 
D2    =   outer    wall    diameter   of    tube 
D;l   =    inner    wall    diameter  of   tube 
k   =   thermal  conductivity   of    wall 
L   =    length    of    heated    surface 


let 


m  Qc    &n    (D2/Pi) 

2     TT     k    L 


(E.3) 
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According    to   Kline      and    HcClintock,       the    uncertainty      in    the 
heat-transfer    coefficient    is: 


6h 
~h" 


6T 


wi 


wo 


-   T 


sat 


+ 


(E.4) 


5F 


wo 


-   T 


sat 


+ 


ST 


sat 


'WO 


-   T 


sat 


1/2 


by  neglecting  the  error  from  the  logarithmic  term,  because 
it  is  small  compared  to  the  other  terms,  the  uncertainty  of 
the    Fourier   term    (F)     can    be    estimated   as: 


and 


6F   «    F 


1/2 


(E.5) 


0C    =    qc    7T    D2    L 


(E.6) 


which   has  an   uncertainty  of 


60, 


1/2 


(E.7) 


Table  5  lists  the  various   terms  of  equations  (E.1)   through 
(E.7)   assuming  all  8  wall  thermocouples  are  used  to  get  the 


average  wall   temperature  f 


wi 


The      data-reduction    program 


(DRP2)       used   this      method      to      calculate      the      heat- transfer 
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TABLE  5 
Uncertainty  Analysis  Terms  Using  8  Thermocouples 


File  WH79         WH79        HF101        HF101 

Heat  Flux     37  kW/ai2     5  xw/m2     37  kfl/m2     5  kW/m2 

6QC 


nc  0.007  0.009  0.006  0.007 

6k 

~k  0.15  0.15  0.33  0.33 


5L 
~L 


0.0005  0.0005  0.0005  0.0005 


F   term     (OC)            0.201                  0.024                1.316  0.170 

SF     (OC)                       0.030                   0.004                0.434  0.056 

0.79  -1.43 

-2.21  -2.18 

0.290  0.  126 

0.259  0.097 

0.015  0.019 

0.008  0.003 

0.389  0.160 

22300  8400 


T„i  <oc> 

1  1.46 

5.63 

Tsat      <°C> 

-2.21 

-2.23 

«T„i 

0.032 

0.073 

"wo    "   Tsat 

6F 

0.002 

0.0005 

T         -    T 
Jwo        1sat 

^sat 

0.003 

0.008 

^wo   "     ^sat 

6qc 

0.003 

0.007 

6h 
~h 

0.032 

0.074 

h    (W/m2    K) 

2660 

540 

I 
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coefficient-  Table  6  lists  the  same  terms  assuming  only  the 
center  4  wall  thermocouples  (2,  3,  5,  and  6)  are  used  to 
calculate  the  heat-transfer  coefficient.  Comparing  Tables  5 
and  6  shows  the  effect  of  the  axial  and  circumferential  wall 
temperature  distributions  on  the  uncertainty  of  the 
heat-transfer  coefficient. 

The  constraining  error  of  the  smooth  tube  is  the  uncer- 
tainty in  the  wail  temperature-  Removing  the  effect  of  the 
axial  wall  temperatue  distribution  reduces  the  uncertainty 
by  1  percent-  The  large  wall  superheat  of  the  smooth  tube 
contributes  to  the  small  magnitude  of  the  uncertainty  terms. 
The  constraining  error  of  the  High  Flux  tube  is  also  the 
uncertainty  in  the  wall  temperature,  but  the  uncertainty  of 
the  thermal  conductivity  (part  of  the  F  term)  is  about  the 
same  magnitude  an  results  in  the  larger  overall  uncertainty 
of  the  High  Flux  data-  The  small  wall  superheats  of  the 
High  Flux  tube  also  amplify  the  magnitudes  of  uncertainty 
terms-  Removing  the  effect  of  the  axial  wall  temperature 
distribution  makes  the  uncertainty  in  the  thermal  conduc- 
tivity of  the  copper-nickel  High  Flux  tube  the  constraining 
uncertainty.  The  axial  temperature  distribution  is  respon- 
sible for  about  25  percent  of  the  wall  uncertainty  term,  but 
again  the  combined  effect  of  the  uncertainty  of  the  thermal 
conductivity,  and  low  wall  superheats,  does  not  make  it 
wholely  responsible  for  the  large  uncertainty  of  the  High 
Flux  data.  More  accurate  data  could  be  obtained  on  the  High 
Flux  surface  by  using  a  solid  copper  tube  without  a  large 
uncertainty  in  the  wall  resistance  or  an  axial  wall 
temperature  variation. 
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TABLE    6 
Uncertainty   Analysis  Teras   Osing   Center    4   Thermocouples 


File  WH79  WH79  HF101  HF101 

Heat    Flux  37    kW/m2  5    kW/m2  37    kW/m2  5    kW/m2 

6f>„ 


Qc 

0.007 

0.009 

0.006 

0.007 

5k 

T 

0.15 

0.  15 

0.33 

0.33 

6L 

~L 

0.0005 

0.0005 

0.0005 

0.0005 

F  term    (°C)  0.201 


0.024 


1.316 


0.  170 


6F    (OQ 


0.030 


0.0  04 


0.434 


0.056 


T         (°C) 

wi 


11.34 


5.80 


0.44 


-1.49 


T  (°C) 

sat 


5T 


wi 


T 

"   Tsat 

6F 

T 
iwo 

Tsat 

^sat 

T 
Lwo 

Tsat 

6qc 

-2.21 


0.015 


0.002 


0.003 


0.003 


0.016 


-2.23 


0.063 


0.0004 


0.008 


0.007 


0.064 


-2.21 


0.213 


0.325 


0.019 


0.008 


0.389 


-2.  18 


0.097 


0.  108 


0.021 


0.003 


0.  146 


h    (W/m2    K)  2690 


530 


27800 


9280 
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